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FOREWORD 


This  report  presents  results  of  work  conducted  by  the  Boeing  Company,  Seattle, 
Washington,  under  Air  Force  Contract  F33615-77-C-3054  "Application  of  the  EASY 
Dynamic  Program  to  the  Analysis  of  Air  Cushion  Systems  on  Aircraft"  during  the 
period  from  15  April  1977  to  1  June  1979.  This  contract  was  conducted  under 
the  sponsorship  of  the  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base,  Ohio  with  Mr.  Peters  Skele  and  Lt.  D.  L.  Fischer  as  project 
engineers. 

This  report  is  comprised  of  three  volumes. 

Volume  I  -  Component  Mathematical  Models 

Volume  II  -  Component  Computer  Programs 

(Parts  I  &  II) 

Volume  III  -  Description  of  Simulations 

In  addition,  a  User's  Manual  (Reference  1)  has  been  written  to  provide  a 
concise  reference  for  day  to  day  usage. 

The  results  presented  were  developed  by  the  Boeing  Aerospace  Company.  The 
program  managers  were  A.  J.  P.  Lloyd,  H.  H.  Straub  and  J.  R.  Kilner.  The 
principal  investigators  were  M.  K.  Wahi,  G.  S.  Duleba,  J.  R.  Kilner  and 
P.  R.  Perkins. 
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SECTION  I 


INTRODUCTION 


The  EASY-ACLS  Model  Generation  and  Analysis  Program  is  documented  in  three 
main  volumes.  Volume  I  contains  the  mathematical  models  for  the  standard 
components  of  the  program  library.  This  includes  derivation  of  the  dynamic 
equations  of  motion  for  aircraft,  trunk-cushion  landing  systems,  an  air  bag 
skid  system  and  an  arresting  gear  system,  and  a  discussion  of  required  inputs 
and  outputs  for  component  models.  Volume  II,  this  volume,  contains  a  detailed 
description  of  all  the  component  subroutines,  including  macro-flowcharts  and 
detailed  listings.  In  addition,  the  standard  subroutines  and  functions  called 
by  the  component  subroutines  are  described,  and  sample  outputs  are  shown. 
Volume  III  contains  a  description  of  the  development,  results  and  conclusions 
of  EASY  computer  simulations  accomplished  during  this  contract.  These  consist 
of  landing,  takeoff,  free  flight  and  drop  test  simulations.  A  summary 
document  (Reference  1)*  has  also  been  prepared  as  a  User's  Manual. 

In  the  EASY  program,  a  user  defines  the  system  to  be  analyzed  by  specifying 
the  individual  components  and  their  interconnections.  Each  of  these 
components  is  a  subroutine  of  the  main  program,  and  each  has  input  data 
requirements  which  must  be  supplied  by  the  user.  In  addition  to  the  library 
of  component  subroutines,  there  is  also  a  set  of  standard  subroutines  and 
functions  which  are  called  automatically  as  they  are  required,  and  for  which 
the  user  need  not  supply  any  input  data.  Such  standard  functions  and 
subroutines  are  described  in  Section  3  of  this  volume-  In  addition,  Section  3 
also  contains  a  description  of  certain  linear  routines  such  as  first  and 
second  order  transfer  functions,  time  delays,  general  controllers,  and 
switches.  Subroutines  associated  with  particular  ACtS  components  are 
described  in  Section  4. 


*References  are  located  in  Volume  II,  Part  II. 
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For  each  function  or  subroutine  described  in  Sections  3  and  4,  the  theory  and 
data  sources  are  given  or  the  appropriate  Section  of  Volume  I  is  referenced. 
Macro- flowcharts  and  some  micro-flow  charts  showing  the  main  decision  points 
and  calculation  procedures  are  shown.  For  all  functions  and  subroutines  a 
listing  is  given.  Included  in  the  listing  is  a  detailed  description  of  the 
input/output  list.  The  inputs  and  outputs  are  divided  into  the  following 
classifications: 

o  Output  State 
o  Output  Derivative 
o  Integrator  Control 
o  Output  Variable 
o  Input  Variable 
o  Input  Parameter 

These  terms  are  discussed  fully  in  Reference  2,  but  will  be  reviewed  briefly 
here  in  order  to  assist  in  understanding  the  subroutines  and  functions.  An 
Output  State  is  a  variable  whose  derivative  (Output  Derivative)  is  calculated 
within  the  subroutine.  The  Integrator  Control  is  an  integer  whose  value  is 
set  by  the  main  program.  The  Output  Derivative  is  calculated  within  the 
subroutine  unless  the  integrator  is  "frozen"  by  setting  the  Integrator  Control 
equal  to  zero.  For  every  Output  State  there  is  one  Output  Derivative  and  one 
Integrator  Control.  Output  Variables,  as  the  name  implies,  are  variables 
which  are  calculated  within  a  function  or  subroutine  and  are  then  output  from 
that  function  or  subroutine.  Similarly,  input  variables  are  variables  which 
die  to  a  function  or  subroutine.  Input  Parameters  are  also  input  to  a 

function  or  subroutine,  but  are  generally  constant  during  a  simulation.  For 
example,  the  inlet  flow  rate  or  temperature  to  a  component  would  be  an  Input 
Variable,  whereas  a  controller  gain  or  duct  length  would  be  an  Input 
Parameter.  In  general,  the  program  user  must  supply  data  inputs  for  all  Input 
Parameters;  however  many  Input  Variables  will  be  automatically  satisfied  by 
the  interconnections  between  adjacent  components.  This  is  discussed  more 
fully  in  Reference  2. 


The  functions  and  subroutines  described  in  Sections  3  and  4  contain  error 
controls  which  are  of  two  types,  namely  for  fatal  and  non-fatal  errors.  The 
latter  type  are  by  far  the  more  common,  and  occur  whan,  for  example,  fluid 
properties  are  sought  outside  the  range  for  which  valid  data  are  available. 

It  is  desirable  that  non-fatal  error  diagnostic  messages  should  not  be  printed 
everytime  they  occur;  for  example,  during  the  course  of  finding  a  steady  state 

solution,  it  is  of  no  significance  if  invalid  data  are  used  provided  the  final 

steady  state  point  is  based  on  valid  data.  The  printout  of  warning 
diagnostics  is  controlled  by  the  parameter  IERR  {=  1  for  printout)  which  is 
set  by  the  main  program  and  is  communicated  through  the  common  block  ERMESS. 
For  errors  which  are  of  sufficient  significance  to  warrant  a  program  stop,  the 

parameter  IFATAL  is  set  equal  to  I  within  the  subroutine  or  function  and  is 

communicated  back  to  the  main  program  through  the  same  common  block  ERMESS. 

Other  common  blocks  which  are  used  in  the  functions  and  subroutines  are: 

CTIME  -  For  time 

CIO  -  For  read  and  write  control 
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SECTION  II 


I 

SUMMARY 


The  EASY  Model  Generation  and  Analysis  program  is  documented  in  three  main 
volumes.  Volume  I  contains  the  mathematical  models.  Volume  II  (this  volume) 
the  component  subroutines,  and  Volume  III  the  description  of  various  dynamic 
simulations.  A  User's  Manual  (Reference  1)  has  also  been  prepared  to  provide 
a  convenient  reference  guide. 

This  volume  contains  a  detailed  description  of  the  component  subroutines,  and 
of  other  standard  functions  and  subroutines  which  are  called  by  the  component 
subroutines.  Flowcharts  and  listings  are  given,  including  a  full  description 
of  the  input/output  lists.  Theory  and  data  sources  are  given,  or  the  relevant 
section  of  Volume  I  is  referenced.  Sample  output  is  included  where 
appropriate.  The  procedures  for  running  single  components  on  the  EASY  program 
are  briefly  reviewed,  and  guidelines  to  assist  users  in  writing  their  own 
subroutines  are  described. 
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SECTION  III 


STANDARD  FUNCTIONS  AND  SUBROUTINES 


3.1  Introduction 

The  EASY  Model  Generation  and  Analysis  program  contains  a  number  of  standard 
functions  and  subroutines  which  have  been  divided  into  six  categories,  namely: 

Property  functions  and  subroutines  -  Section  3.2 

Flow  functions  and  subroutines  -  Section  3.3 

Transfer  functions  and  subroutines  -  Section  3.4 

Miscellaneous  subroutines  -  Section  3.5 

Table  look  up  functions  -  Section  3.6 

Initial  condition  functions  and  subroutines  -  Section  3.7 

Supporting  subroutines  for  Foster-Miller  -  Section  3.8 

trunk  component 


The  property,  the  flow,  the  table  look  up,  and  the  inital  condition  functions 
and  subroutines  are  all  called  by  component  subroutines,  and  their  input  data 
requirements  are  automatically  satisfied  by  the  calling  subroutine.  The 
transfer  functions  are  generally  used  to  model  part  or  all  of  a  component,  and 
are  therefore  included  in  the  library  of  standard  components  (see  Reference  1). 
However,  as  they  are  general  dynamic  models  and  are  not  associated 
particularly  with  any  one  type  of  ACLS  component,  they  are  included  in  this 
section  rather  than  in  the  component  subroutines  discussed  in  Section  4  of 
this  volume.  The  miscellaneous  functions  include  routines  for  providing  input 
data  to  the  models  in  tabular  form,  and  for  changing  variable  or  parameter 
inputs  during  a  simulation.  Table  1  lists  all  the  standard  functions  and 
subroutines  and  briefly  describes  their  purpose.  A  full  description  of  each 
routine  is  given  in  the  following  sections. 
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Table  1  STANDARD  FUNCTIONS  AND  SUBROUTINES 


FUNCTION  OR 


SECTION 

NO. 

SUBROUTINE 

NAME 

PURPOSE 

3.2.1 

SHCP 

Specific  heat  of  humid  air 

3.2.2 

PROP 

Property  of  fluids 

3.2.3 

PRND 

Prandtl  number  of  fluids 

3.3.1 

FNFLOW 

Flow  rate  in  a  duct  given  pressure  and 
loss  factor,  using  tabular  data 

3.3.2 

FSFLOW 

Flow  rate  in  a  duct  given  pressure 
and  loss  factor,  using  tabular  data 

3.3.3 

PERF 

Trunk  element  orifice  areas  for  flow 
to  cushion  and  atmosphere 

3.3.4 

PERFB 

Air  bag  element  orifice  areas  for 
flow  to  atmosphere 

3.3.5 

AMACH 

Mach  number  in  a  duct 

3.3.6 

RENVX 

Reynolds  number 

3.3.7 

HI 

Convective  heat  transfer  coefficients 
for  flow  in  ducts  and  across  cylinder 

3.3.8 

VI X 

Flow  across  butterfly,  gate  and  globe 
val ves 

3.4.1 

LA 

First  order  lag  (time  constant  form) 

3.4.2 

LG 

First  order  lag  (pole  form) 

3.4.3 

LL 

Lead  lag  (time  constant  form) 

3.4.4 

LE 

Lead  lag  (pole-zero  form) 

3.4.5 

TF 

Second  order  transfer  function  with  first 
order  numerator 
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Table  1  STANDARD  FUNCTIONS  AND  SUBROUTINES 


FUNCTION  OR 


SECTION 

NO. 

SUBROUTINE 

NAME 

PURPOSE 

3.4.6 

TZ 

Second  order  transfer  function  with  second 
order  numerator 

3.4.7 

RG 

Rate  Gyro  Dynamics  and  Saturation 

3.4.8 

FG 

F  ight-ground  controller 

3.4.9 

AP 

Autopilot  pitch  controller  (Jindivik) 

3.4.10 

AR 

Autopilot  roll  controller  (Jindivik) 

3.4.11 

IT 

Integrator  with  saturation 

3.5.1 

SW 

Switch  control  for  one  variable 

3.5.2 

sx 

Switch  control  for  two  variables 

3.5.3 

SY 

Switch  control  for  three  variables 

3.5.4 

SZ 

Switch- control  for  four  variables 

3.5.5 

TA 

Tauular  input  routine  for  two 
variables 

3.5.6 

TB 

Tabular  input  routine  for  four  variables 

3.5.7 

AF 

Analytical  function  generator 

3.5.8 

MA 

Multiply  and  add 

3.5.9 

FU 

Arbitrary  function  generator 

3.5.10 

MC 

Multiply  and  add  with  three  variables 

3.5.11 

SA 

Saturation 

3.5.12 

SB 

Saturation  with  dead  band 

3.5.13 

MB 

Multiply,  add  and  divide 

3.5.14 

FV 

Two  dimension  function  generator 
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Table  1  STANDARD  FUNCTIONS  AND  SUBROUTINES 


FUNCTION  OR 


SECTION 

NO. 

SUBROUTINE 

NAME 

PURPOSE 

3.5.15 

RA 

Random  number  generator  function 

3.5.16 

RN 

Random  number  generator 

3.5.17 

S2 

Sum  forces  and  moments  (2  sets  of 
inputs) 

3.5.18 

S3 

Sum  forces  and  moments  (3  sets  of 
incuts) 

3.5.19 

TG 

Transform  engine  thrust  into  body  axis 

3.5.20 

TR 

Transform  vectors  body  to  earth  axis 

3.5.21 

X? 

Transform  angular  rates 

3.5.22 

XT 

Transform  torques 

3.5.1 

T3L1 

Linear  interpolation  (one  ind. 
variable) 

3.6.2 

73L2 

Linear  interpolation  (two  ind. 
variables) 

3.6.3 

TBLU1 

Polynomial  interpolation  (one  ind. 
variable) 

3.5.4 

TBLU2 

Polynomial  interpolation  (two  ind. 
variables) 

3.6.5 

TBLU3 

Polynomial  interoolation  (three  ind. 
variables) 

3.6.5 

ET32 

Linear  interpolation  (two  independent  variables) 
Supports  TS 

^  •  D  •  / 

EIB3 

Linear  interoolation  (three  independent  variables 
Supports  TS 

7  » 

O  .  ..  i 

*  ^ 

i  L 

Initial  condition  oackage  (trunk' 
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Table  1  STANDARD  FUNCTIONS  AND  SUBROUTINES 


FUNCTION  OR 


SECTION 

NO. 

SUBROUTINE 

NAME 

PURPOSE 

3.7.2 

ICB 

Inital  condition  package  (air  bag) 

3.7.3 

ICFS 

Free  shape  equations  (trunk) 

3.7.4 

ICFSB 

Free  shape  equations  (air  bag) 

3.7.5 

I  CIS 

Loaded  shape  equations  (trunk) 

3.7.6 

ICLSB 

Loaded  shape  equations  (air  bag) 

3.7.7 

KINK 

Kink  wave  angle  equations  (arresting 
cable) 

3.7.8 

RES 

Arresting  cable  strain  equations 

3.7.9 

TERRA 

Terrain  model 

3.7.10 

VPRIN8 

Air  bag  arrays  print  control 

3.7.11 

VPRINT 

Trunk  arrays  print  control 

3.7.12 

ELAS 

Initial  condition  package  for  elastic  trunk  TS 

3.7.13 

ELFX 

Incomplete  elliptic  integral  of  the  2nd  kind 

3.7.14 

ELKX 

Complete  elliptic  integral  of  the  2nd  kind 

3.7.15 

ELWR 

Prints  properties  of  elastic  trunk  TS 

3.7.16 

ENDFS 

Equations  for  elastic  trunk-end  element-free  shape 

3.7.17 

ENDLS 

Equations  for  elastic  trunk-end  element-loaded 
shape 

0  7  10 

*J  •  /  m  iU 

SIDEFS 

Equations  for  elastic  trunk-side  element-free  shape 

3.7,19 

SIDELS 

Equations  for  elastic  trunk-side  element-loaded 
shape 

3.7.20 

XXPRT 

Prints  elastic  trunk  parameters  during  simulation 

If 

& 

sp 

y*yS 

IS 
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Table  1  STANDARD  FUNCTIONS  AND  SUBROUTINES 


FUNCTION  OR 


SECTION 

NO. 

SUBROUTINE 

NAME 

PURPOSE 

3.8.1 

CDVCHP 

Updates  cushion  volume-pressure  parameters  (FM)* 

3.8.2 

CLRNCE 

Ground  clearance  for  trunk  segments  (FM)* 

3.8.3 

COQRDN 

Ground  coordinates  of  trunk  segments  (FM)* 

3.8.4 

DYNFAN 

Fan  dynamic  characteristics  (FM)* 

3.8.5 

FLOW 

Dynamic  flow  in  ACLS  (FM)* 

3.8.6 

FMFAN 

Fan  static  characteristics  (FM)* 

3.8.7 

FMWRIT 

Prints  ACLS  input  data  (FM)* 

3.8.8 

FORCE 

ACLS  forces  (FM)* 

3.8.9 

HYCURV 

Side  trunk  height  (FM)* 

3.8.10 

OUTFM 

Updates  variables  during  simulation  (FM)* 

3.8.11 

PARAMS 

Def vault  values  for  ACLS  parameters  (FM)* 

3.8.12 

PROFILE 

Defines  ground  elevation  (FM)* 

3.8.13 

ROTATE 

Vehicle-to-initial  vector  transformation  (FM)* 

3.8..  14 

SEGMNT 

Divides  trunk  into  segments  (FM)* 

3.8.15 

SHAPE 1 

Initial  estimates  of  ACLS  areas  and  volumes  (FM) 

3.8.16 

SHAPE 2 

ACLS  areas  and  volumes  (FM)* 

3.3.17 

STATIC 

Load  maps  and  static  equilibrium  conditions  (FM) 

3.8.18 

STEQU 

Defines  ACLS  state  derivatives  (FM)* 

3.8.19 

TRUNK 

Trunk  cross  section  parameters  (FM)* 

3.8.20 

VALVE 

Relief  valve  orifice  area  (FM)* 

*FOsTER-MILLER  ACLS  MODEL 
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3.2  Property  Functions  and  Subroutines 


3.2.1  Function  SHCP 

The  purpose  of  function  SHCP  is  to  calculate  the  specific  heat  of  an  air-water 
vapor  mixture.  The  function  is  called  by  many  of  the  component  subroutines 
for  calculating  specific  heat,  and  hence  the  gas  constant  and  the  ratio  of 
specific  heats.  The  method  is  described  in  Reference  2,  Volume  I,  Appendix  A 
and  uses  property  data  from  References  3  and  4.  The  function  is  valid  for 
temperatures  of  300-1600°R,  and  a  warning  diagnostic  message  is  printed  if  the 
temperature  is  outside  this  range.  Table  88  shows  a  flowchart  and  Table  214 
gives  a  detailed  listing.  Required  inputs  are  the  temperature  (T)  and 
specific  humidity  (SH). 

3.2.2  Function  PROP 

The  purpose  of  function  PROP  is  to  calculate  the  properties  of  several  fluids. 
The  properties,  all  of  which  are  assumed  to  be  functions  of  temperature  only, 
are  specific  heat,  viscosity,  thermal  conductivity  and  density.  The  fluids 
included  in  the  routine  are: 

o  Dry  air 

o  Water  at  saturation  pressure 
o  60%/40%  Ethylene  Glycol /water 
o  Heat  transport  fluid  FC-75 
o  Fuel  JP-4  (MIL-F-5624) 
o  Heat  transport  fluid  DC-331 
o  Hydraulic  fluid  (MIL-H-83282) 
o  Hydraulic  fluid  (MIL -H-5606 ) 

The  function  uses  polynomial  approximations  to  property  data  obtained  from 
References  3,  5  and  6.  Table  73  shows  a  flowchart  and  Table  199  gives  a 
detailed  listing.  Limitations  on  the  allowable  temoerature  ranoe  fnr  each 
fluid  are  shown  m  the  listing.  A  warning  diagnostic  is  printed  if 
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temperatures  are  outside  the  allowable  range.  The  density  calculation  for  dry 
air  is  invalid,  and  the  perfect  gas  equation  should  be  used  instead.  Required 
inputs  are  integers  identifying  the  fluid  and  the  required  property,  and  the 
absolute  temperature. 

3.2.3  Function  PRND 

The  purpose  of  function  PRND  is  to  calculate  the  Prandtl  number  of  fluids. 

The  function  calls  PROP  (Section  3.2.2)  to  calculate  the  specific  heat, 
thermal  conductivity  and  viscosity,  with  the  exception  that  the  specific  heat 
of  moist  air  is  obtained  from  SHOP.  As  discussed  in  Reference  2,  Volume  I, 
Appendix  A  values  of  thermal  conductivity  and  viscosity  for  dry  air  can  be 
used  in  the  calculation  of  Prandtl  number  of  an  air/water  vapor  mixture  with  a 
neglible  effect  on  accuracy.  No  diagnostic  messages  are  included  in  PRND,  but 
will  be  called  from  PROP  and  SHOP  (Section  3.2.1)  if  the  temperature  lies 
outside  the  valid  range  for  the  particular  fluid.  Table  71  shows  a  flowchart 
and  Table  197  gives  a  detailed  listing.  Input  requirements  are  an  integer 
identifying  the  fluid  (as  in  PROP),  the  absolute  temperature  and  the  specific 
humidity.  The  latter  input  is  used  for  the  case  of  air  only. 

3.3  Flow  Functions  and  Subroutines 
3.3.1  Subroutine  FNFLOW 

The  purpose  of  subroutine  FNFLOW  is  to  calculate  the  flow  rate  and  Chester 
Smith  function  for  flow  in  a  constant  area  section,  given  the-  upstream  and 
downstream  pressures,  temperature,  effective  area  and  K  (pressure  loss) 
factor.  The  method  is  described  in  Reference  2,  Volume  I,  Appendix  F.  The 
Chester  Smith  function  is  stored  in  a  data  array  as  a  function  of  pressure 
ratio,  .'he  subroutine  can  also  be  used  to  calculate  the  flow  across  an 
orifice,  in  which  case  the  K  factor  should  be  set  equal  to  one.  The  input 
value  of  K  -'actor  must  a'  .’ays  be  greater  than  zero.  In  the  event  that  the 
input  value  of  downstream  pressure  exceeds  the  upstream  pressure,  the 
subroutine  ca'culates  negative  flow  rates  and  Chester  Smith  functions. 
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A  flowchart  for  the  subroutine  is  shown  as  Table  37  and  Table  163  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.3.2  Subroutine  FSFLOW 

The  purpose  of  subroutine  FSFLOW  is  the  same  as  that  of  FNFLOW  except  an 
additional  output  parameter  "SFN"  has  been  added.  Parameter  SFN  is  the  slope 
of  compressible  flow  factor  FN  with  respect  to  the  pressure  ratio  P1/P2. 
Subroutine  FSFLOW  is  used  exclusively  for  calculating  cushion  to  atmosphere 
flow  rate  in  the  trunk  subroutine  TK. 

A  flowchart  for  the  subroutine  is  shown  as  Table  41  and  Table  167  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.3.3  Subroutine  PERF 

This  is  a  special  purpose  subroutine  called  only  by  the  trunk  component  TK. 

The  purpose  of  subroutine  PERF  is  to  calculate  orifice  areas  for  flow  to 
atmosphere  and  to  cushion  through  the  lubrication  holes  provided  on  the 
underside  of  a  trunk.  Thirteen  combinations  of  performation  arrangements  are 
treated  and  are  functions  of  free  or  loaded  trunk  shape  and  perforations  being 
on  the  flattened  or  free  parts  of  the  trunk. 

A  flowchart  for  the  subroutine  is  shown  as  Table  69  and  Table  195  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.3.4  Subroutine  PERFB 

This  is  also  a  special  purpose  subroutine  called  only  the  the  Air  Bag 
Component  AB.  Its  purpose  is  same  as  that  of  the  subroutine  PERF  except  the 
perforation  arrangements  which  are  limited  to  six  and  there  is  no  cushion 
pressure. 
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A  flowchart  for  the  subroutine  is  shown  as  Table  70  and  Table  196  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.3.5  Function  AMACH 

The  purpose  of  function  AMACH  is  to  calculate  the  Mach  number  of  moist  air 
flowing  in  an  element  of  known  cross-sectional  area,  given  the  pressure, 
temperature,  flow  rate  and  specific  humidity.  The  method  is  described  in 
Reference  2,  Volume  I,  Appendix  F.  The  function  f(M)  defined  in  equation  F-2 
is  stored  in  two  data  arrays.  The  first  array  covers  a  range  of  values  of 
f (M)  from  zero  through  0.50  at  increments  of  .02,  corresponding  to  Mach 
numbers  of  zero  through  .632.  The  second  array  covers  the  range  .50  through 
.58  at  increments  of  .01,  corresponding  to  Mach  numbers  of  .632  through  1.0. 

A  warning  diagnostic  message  is  printed  if  the  calculated  value  of  f(M)  lies 
outside  the  range  covered  by  the  data  arrays. 

A  flowchart  for  the  function  is  shown  as  Table  5  and  table  131  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.3.6  Function  RFNVX 

The  purpose  of  function  RENVX  is  to  calculate  the  Reynolds  number  of  fluid 
flowing  in  a  duct  of  circular  cross-section,  given  the  flow  rate,  temperature 
and  duct  diameter.  The  function  calls  function  PROP  (Section  3.2.2)  to 
determine  the  fluid  viscosity  at  the  given  temperature.  As  discussed  in 
Reference  2,  Volume  I,  Appendix  A,  values  of  viscosity  for  dry  air  can  be  used 

for  the  calculation  of  Reynolds  number  of  an  air/water  vapor  mixture  with 

negligible  effect  on  accuracy.  The  function  output  is  the  absolute  value  of 

Reynolds  number,  irrespective  of  the  sign  of  the  flow  rate. 

A  flowchart  for  the  function  is  shown  as  Table  76  and  Table  202  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 
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3.3. 7  Function  HI 


The  purpose  of  function  HI  is  to  calculate  the  convective  heat  transfer 
coefficient  for  fluid  flow  in  ducts  and  for  flow  across  cylinders.  The  latter 
capability  is  primarily  of  significance  for  determining  the  response  rate  of 
temperature  sensors.  The  function  uses  conventional  correlations  for  Nusselt 
number  against  Reynolds  rumber  and  Prandtl  number  for  transitional  and 
turbulent  flow  in  ducts,  and  for  Nusselt  number  against  Graetz  number  f'  . 
laminar  flow  in  ducts.  For  flow  across  cylinders,  a  correlation  of  N'^selt 
number  against  Reynolds  number  is  used  for  air,  and  Nusselt  number  against 
Reynolds  number  and  Prandtl  number  for  liquids.  The  complete  equations  and 
references  are  given  in  Reference  2,  Volume  I,  Appendix  F.  Fluid  properties 
are  obtained  from  function  PROP  (Section  3.2.2). 

A  flowchart  for  the  function  is  shown  as  Table  46  and  Table  172  gives  a 
detailed  listing  which  includes  the  full  input/output  list.  Figure  1  shows 
heat  transfer  coefficients  for  air  flow  in  ducts  and  across  cylinders. 

3.3.8  Subroutine  VLX 

The  purpose  of  subroutine  VLX  is  to  calculate  the  weight  flow  of  air  across 
butterfly,  gate  and  globe  type  valves.  For  butterfly  and  globe  valves,  the 
geometric  flow  area  is  first  calculated,  and  the  effective  area  is  then 
determined  assuming  discharge  coefficients  of  0. &/  and  0.80,  respectively. 

For  gate  valves  the  flow  is  calculated  using  a  table  of  K  factor  versus 
fractional  opening.  A  full  description  of  the  method  and  equations  is  given 
in  Reference  2,  Volume  I,  Section  3,7.  In  the  event  that  the  downstream 
pressure  is  greater  than  the  upstream  value,  the  calculated  value  of  flow  will 
be  negative.  For  globe  valves,  the  input  value  of  poppet  diameter  must  be 
greater  than  or  equal  to  the  poppet  seat  diameter. 

The  valve  opening  (degrees  or  fractional  opening)  is  an  input  to  subrovfine 
VLX.  In  a  system  analysis,  the  valve  opening  will  generally  be  an  output 
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variable  from,  or  calculated  within,  a  control  component  such  as  a  pressure 
regulator  or  temperature  controller.  During  the  course  of  searching  for  a 
steady  state  solution,  it  is  possible  that  the  valve  opening  will  lie  outside 
the  normal  range.  For  example,  a  butterfly  valve  opening  of  -10°  or  105° 
might  be  input  to  VLX.  Although  such  a  value  is  clearly  not  a  valid  in  a 
steady  state  solution,  it  is  desirable  that  subroutine  Vi.X  calculate 
reasonable  values  of  flow  rate  at  such  values  in  order  to  permit  convergence 
to  the  correct  steady  state  solution.  The  method  incorporated  into  VLX  is 
illustrated  in  Figure  2(a)  for  butterfly  valves,  and  Figure  2(b)  for  globe  and 
gate  valves. 

A  flowchart  for  the  subroutine  is  shown  as  Tabic  120  and  Table  246  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.4  Transfer  Functions  and  Subroutines 

3,4.1  Subroutine  LA 

The  purpose  of  subroutine  LA  is  to  simulate  a  first  order  lag  defined  as: 

F0  _  Gain  3.4-1 

FIN  (1  +  TC  x  s) 

where 

FO  is  the  output  GAIN 

FIN  is  the  input  TC 

and  s  is  the  Laplace  operator 

The  time  constant  TC  must  not  equal  zero. 

Table  56  shows  a  flowchart  for  the  subroutine  and  Table  182  gives  a  detailed 
listing  which  includes  the  full  input/output  list.  Figure  3(a)  shows  a 
typical  calculated  response  to  a  unit  step  change  in  input,  which  agrees  with 
the  analytical  response  given  by 


is  the  gain 

is  the  time  constant 
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AMPLITUDE 


FO  =  (1  -  e"t/TC)  x  GAIN  3.4-2 

where  t  is  time.  Figure  3(b)  shows  the  frequency  response. 

3.4.2  Subroutine  LG 

The  purpose  of  subroutine  LG  is  to  simulate  a  first  order  lag  defined  by: 


where 


FO  ZO 

m  =  im 

FO  is  the  output 
FIN  is  the  input 


ZO  is  the  zero  location 
PO  is  the  pole  location 


3.4.3 


and  s  is  the  Laplace  operator. 

A  flowchart  for  the  subroutine  is  shown  as  Table  58  and  Table  184  gives  a 
detailed  listing. 

3.4.3  Subroutine  LL 

The  purpose  of  subroutine  LL  is  to  simulate  a  lead-lag  defined  by: 


FO 

FIN 


where 


GAIN 


1  +  TCI  x  s) 
1  +  TC2  x  s) 


3.4-4 


FO  is  the  output  GAIN  is  the  gain 

FIN  is  the  input  TCI  and  TC2  are  time  constants 


and  s  is  the  Laplace  operator. 

The  time  constant  TC2  must  not  equal  zero. 


-3* 
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Table  59  shows  a  flowchart  for  the  subroutine  and  Table  185  gives  a  detailed 
listing  which  includes  the  full  input/output  list.  Table  4(a)  shows  a  typical 
calculated  response  to  a  unit  step  change  in  input,  which  agrees  with  the 
analytical  response  given  by 


FO  „  .  (TCI  -  TC2)  -t/TC2*  „  r,TU 

m  =  vi  + — K2 —  e  }  x  GAIN 

where  t  is  time.  Figure  4(b)  shows  the  frequency  response. 


i.4-5 


3.4.4  Subroutine  LE 

The  purpose  of  subroutine  LE  is  to  simulate  a  first  order  lead- lag  defined  by: 


FO 

rm 


GAIN 


3.4-6 


where 


FO 

is 

the  output 

GAIN 

is 

FIN 

is 

the  input 

ZO 

is 

s 

is 

the  Laplace 

operator.PO 

is 

A  flowchart  for  the  subroutine  is  shown  as  Table  57  and  Table  183  gives  a 
detailed  listing. 

3.4.5  Subroutine  TF 

The  purpose  of  subroutine  TF  is  to  simulate  a  second  order  transfer  function 
with  first  order  numerator  defined  by: 


FO 

Hn 


Z1  x  s  +  ZO 


3.4-7 


S*  +  PI  x  s  +  PO 
where  FO  is  the  outout 
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FIN  is  the  input 

ZP,  Zl,  PO,  and  PI  are  coefficients 
and  s  is  the  Laplace  operator 


Table  109  shows  a  flowchart  for  the  subroutine  and  Table  235  gives  a  detailed 
listing  which  includes  the  full  input/output  list.  Figure  5(a)  shows  the 
calculated  response  to  a  unit  step  change  in  input,  which  agrees  with  the 
analytical  response  given  by 


FO  = 


_Yt  ,B  _YAv  Jr 

Ae  Cos  u>t  +  (  •  ^  -  e  Sin  tot  -  A 


3.4-8 


where 


ZO 

PIT 


B  =  Zl  - 


ZO  x  PI 


ur  *  PO¬ 


TT 

(Pl)2 


PI 

T 


and  t  is  the  time 

Figure  5(b)  shows  the  frequency  response. 


.6  Subroutine  TZ 


The  purpose  of  subroutine  TZ  is  to  simulate  a  second  order  transfer  function 
with  second  order  numerator  defined  by 


FO 

FTTT 


Z2  x  s^  +  Zl  x  s  +  ZO 


3.4-9 


s  +  PI  x  s  +  PO 
where  FO  is  the  output 
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FIN  is  the  input 

ZO,  Zl,  Z2,  PO,  and  PI  are  coefficients 
and  s  is  the  Laplace  operator. 

Table  117  shows  a  flowchart  for  the  subroutine  and  Table  243  gives  a  detailed 
listing  which  includes  the  full  input/output  list.  Figure  6(a)  shows  the 
calculated  response  to  a  unit  step  change  in  input,  which  agrees  with  the 
analytical  response  given  by 
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Figure  6 


Typical  Response  for  Subroutine  TZ  ' 
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A  flowchart  for  the  subroutine  is  shown  as  Table  78  and  Table  204  gives  a 
detailed  listing  which  includes  the  input/output  list. 


3*4.8  Subroutine  FG 


.gsgSag 


nd  =  .28  (e-  ed)  +  .533  (  8-ed) 

where  nd  =  elevator  angle 
§  =  pitch  rate 
Qd  =  pitch  reference  rate 
8  =. pitch  angle 
0d  =  pitch  reference  angle 
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The  elevator  se; vo  transfer  function  is 


JX,  = _ 1  +  .  32S 

nd  (1+.4S) (1+.11S+-0063  s2) 


3.4-12 


where  s  is  the  Laplace  transform  operator. 

A  flowchart  for  the  subroutine  is  shown  in  Table  6  and  Table  132  gives  a 
detailed  listing. 

3.4.10  Subroutine  AR 

The  purpose  of  subroutine  AR  is  to  simulate  JINDIVIK  autopilot  for  roll 
control  and  is  only  a  special  purpose  subroutine. 

The  lateral  autopilot  function  is 


Cd  =  0.196  (*-*d)  +  .42  (*-#d)  +  -2R 
+  Kr  /Rdt  +  .0082jQfRdt 


3.4-13 


where  Cd  =  aileron  angle  contnanded 
$  =  rol  1  rate 

<$d  =  roll  reference  rate  (lOdeg/sec) 
0=  roll  angle 

#d  =  roll  reference  angle 
R  =  yaw  rate 

Kr  =  0.35  for  flight  and  approach 

K  =  0  after  touchdown  (t=0) 
r 
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Cd  (1+.4S) (1+.11S+.0063  s2) 
where  s  is  the  Laplace  transform  operator. 

A  flowchart  for  the  subroutine  is  shown  in  Table  7  and  Table  133  gives  a 
detailed  listing. 

3.4.11  Subroutine  IT 

The  purpose  of  subroutine  IT  is  to  simulate  an  integrator  with  limits  on  the 
upper  and  lower  limits  of  the  output.  Table  54  shows  a  flowchart  of  the 
subroutine  and  Table  180  gives  a  detailed  listing  which  includes  the  full 
input/output  list.  The  values  of  AMA,  AMI  and  GKL  can  be  selected  using  the 
same  logic  described  in  Reference  2,  Volume  II,  Section  3.4-10. 

3,5  Miscellaneous  Subroutines 

3.5.1  Subroutine  SW 

The  purpose  of  subroutine  SW  is  to  provide  the  capability  to  switch  the  value 
of  a  variable  during  the  course  of  a  simulation.  The  subroutine  output  (V01) 
is  initially  set  equal  to  one  or  other  of  the  input  variables  (VA1  and  V81), 
depending  on  the  initial  value  of  the  switch  control  (SW1).  At  time  TCI  the 
output  variable  is  switched,  and  at  time  TC2  (with  TC2  greater  than  TCI)  the 
output  variable  is  switched  back. 

The  subroutine  may  be  used  to  provide  forcing  function  inputs  such  as  step 
changes  or  impulses,  to  the  ooundary  conditions  of  a  system,  as  shown  in 
Figure  7(a)  and  7(b).  ( Arbitrary  time  dependent  inputs  can  be  obtained  from 

subroutines  TA  and  T8  described  in  Sections  3.5.5  and  3.5.6  respectively). 

The  subroutine  can  also  be  used  to  change  the  value  of  an  input  parameter  of  a 
component.  For  example,  the  K  factor  of  a  duct  could  be  changed  during  a 
simulation  in  order  to  determine  the  effect  of  a  partial  flow  blockage;  this 
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Application  of  Switch  Subroutine  SW 


is  illustrated  in  Figure  7(c).  Also  the  subroutine  can  be  used  to  switch  the 
output  signals  from  two  seperate  components  as  inputs  into  a  single  component, 
(Figure  7(d)). 

A  flowchart  for  the  subroutine  is  shown  as  Table  94  and  Table  220  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.5.2  Subroutine  SX 

The  purpose  of  subroutine  SX  is  to  provide  the  capability  to  switch  the  values 
of  two  variables  during  the  course  of  a  simulation.  The  subroutine  is  very 
similar  to  SW  (described  in  Section  3.5.1)  except  that  there  are  two  outputs 
(V01  and  V02),  with  two  inputs  for  each  output  (VA1  and  VB1  for  V01,  VA2  and 
VC2  for  V02).  The  outputs  are  both  switched  at  time  TCI  and  then  switched 
back  at  time  TC2  (with  TC2  greater  than  TCI).  Possible  applications  of 
subroutine  SX  are  as  discussed  in  Section  3.5.1  and  illustrated  in 
Figure  7. 

The  flowchart  for  the  subroutine  is  shown  as  Table  95  and  Table  221  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.5.3  Subroutine  SY 

The  purpose  of  subroutine  SY  is  to  provide  the  capability  to  switch  the  values 
of  three  variables  during  the  course  of  a  simulation.  The  subroutine  is  very 
similar  to  SW  (described  in  Section  3.5.1),  except  that  there  are  three 
outputs  (VOX,  V02,  and  V03)  with  two  inputs  for  each  output  (VA1  and  VB1  for 
V01,  VAZ  and  VB2  for  V02,  and  VA3  and  VB3  for  V03).  The  outputs  are  all 
switched  at  time  TCI  and  then  switched  back  at  time  TC2  (with  TC2  greater  than 
TCI).  Possible  applications  of  subroutine  SY  are  discussed  in  Section  3.5.1 
and  illustrated  in  Figure  7. 

The  flowchart  for  the  subroutine  is  shown  as  Table  96  Table  222  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 
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3.5.4  Subroutine  SZ 


The  purpose  of  subroutine  SZ  is  to  provide  the  capability  to  switch  the  values 
of  four  variables  during  the  course  of  a  simulation.  The  subroutine  is  very 
similar  to  the  SX  (as  described  in  Section  3.5.1)  except  that  there  are  four 
Outputs  (VC1,  V02,  V03,  and  V04)  with  two  inputs  for  each  output  (VA1  and  VB1 
for  V01,  VA2  and  V82  for  V02,  VA3  and  VB3  for  V03,  and  VA4  and  VB4  for  ¥04). 
The  outputs  are  all  switched  at  time  TCI  and  then  switched  back  at  time  TC2 
(with  TC2  greater  than  TCI).  Possible  applications  of  subroutine  SZ  are  as 
discussed  in  Section  3.5.1  and  illustrated  in  Figure  7. 

The  flowchart  for  the  subroutine  is  shown  as  Table  97.  Table  223  gives  a 
detailed  listing  which  includes  the  full  input/output  list. 

3.5.5  Subroutine  TA 

The  purpose  of  subroutine  TA  is  to  provide  the  capability  to  input  up  to  four 
variables  as  functions  of  time.  Inputs  to  the  subroutine  are  four 
one-dimensional  tables  (A2TAB,  B2TAB,  C2TA8,  and  D2TAB)  with  time  as  the 
independent  variable.  The  corresponding  outputs  are  the  variables  A2,  B2,  C2, 
and  02. 

The  value  of  time  is  obtained  within  the  subroutine  from  the  common  block 
CTIME. 

A  typical  application  of  the  subroutine  would  be  to  calculate  engine  bleed  and 
fan  air  pressures  and  temperature  during  a  throttle  burst  or  chop. 
Alternatively  it  can  be  used  to  input  other  time  dependent  variables  such  as 
electrical  heat  loads.  Step  changes  or  impulses  can  be  input  more  efficiently 
using  the  "switch"  subroutines  SW,  SX,  SY  and  SZ  (Sections  3.5.1  through 
3.5.4).  Subroutine  TA  should  be  used  where  there  are  four  or  three  time 
dependent  inputs  (in  the  latter  instance  one  table  is  a  dummy);  for  two  or  one 
time  dependent  inputs,  subroutine  TB  (Section  3.5.6)  is  more  appropriate. 
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A  flowchart  for  the  subroutine  is  shown  as  Table  100  and  Table  226  gives  a 
detailed  listing  which  includes  the  full  input /output  list. 

3.5.6  Subroutine  TB 

The  purpose  of  subroutine  TB  is  to  input  up  to  two  variables  as  a  function  of 
time.  Inputs  to  the  subroutine  are  two  one-dimensional  tables  (A2TAB  and 
B2TAB)  with  time  as  the  independent  variable.  The  corresponding  outputs  are 
the  variable  A2  and  B2,  The  value  of  time  is  obtained  within  the  subroutine 
from  the  common  block  CTIME. 

This  subroutine  is  very  similar  to  TA  (Section  3.5.5)  except  that  it  is 
limited  to  two  input  tables  and  corresponding  outputs,  whereas  TA  has  the 
capability  of  up  to  four  input  tables.  Typical  applications  of  both  TA  and  TB 
are  discussed  in  Section  3.5.5. 

The  flowchart  for  the  subroutine  is  shown  as  Table  101.  Table  227  gives  a 
detailed  listing  for  subroutine  TB,  including  the  full  input/output  list. 

3.5.7  Subroutine  AF 

Tne  purpose  of  subroutine  AF  is  to  generate  several  analytical  functions,  with 
time  as  the  indepedent  variable.  The  output  from  the  subroutine  can  be  used 
as  an  input  to  a  system  in  order  to  determine  the  system  response  to  the 
particular  form  of  excitation.  The  analytical  function  is  selected  by  an 
input  code  as  follows: 
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COD  =  OUTPUT  = 

1  C,  +  C2  Sin  (C3  t  +  C4). 

2  C1  +  C2  Cos  ^C3  t  +  c*)* 

3  C  +  C„  e’V  Sin  (C„t  +  C,). 

1  C  O  4 

4  C1  +  C2  e'C5t  C0S  ^C3t  +  C4'* 


where  C. ,  C0,  C.,,  C,,  and  C,  are  constants  and  t  is  time. 

1  2’  3’  4*  d 

Table  4  shows  a  flowchart  for  the  subroutine,  and  Table  130  gives  a  detailed 
listing  which  includes  the  full  input/output  list. 

3.5.8  Subroutine  MA 

The  purpose  of  Subroutine  MA  is  to  reset  the  value  of  a  variable  according  to 
the  equation 

FO  =  Cl  x  FIN  +  C2  3.5-2 


where 

FQ 

is  the  output 

FIN 

is  the  input 

and 

Cl,  C2 

are  input  constants 

A  typical  application  for  subroutine  MA  would  be  to  simulate  a  symmetrical 
flow  merge,  for  example  where  the  flow  from  two  engine  bleed  systems  merge  to 
flow  into  a  single  air  conditioning  pack.  FGr  such  a  case  FIN  would  be  the 
flow  from  a  single  engine,  FQ  would  be  the  pack  flow.  Cl  would  be  set  equal  to 
2  and  C2  to  zero. 
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Table  60  shows  a  flowchart  for  the  subroutine  and  Table  186  gives  a  detailed 
listing  which  includes  the  full  input/cutput  list. 

3.5.9  Subroutine  FU 

The  purpose  of  subroutine  FU  is  to  provide  the  capability  of  simulating  an 
arbitrary  one-dimensional  table  y  =  f(X).  The  user  has  the  capability  of 
specifying  the  degree  of  interpolation  and  whether  or  not  extrapolation  should 
be  permitted. 

Table  43  shows  a  flo'.rchart  for  the  subroutine  and  Table  169  gives  a  detailed 
listing  which  includes  the  full  input/output  list. 

3.5.10  Subroutine  MC 

The  purpose  of  subroutine  MC  is  to  reset  the  value  of  a  variable  according  to 
the  equations 

FO  =  Cl  x  FIN  +  C2  x  FIO  +  C3  x  FIP  +  C4  3.5-3 

where  FO  is  the  output 

FIN,  FIO,  FIP  are  inlet  variables 
and  Cl,  C2,  C3  and  C4  are  input  constants. 

A  typical  application  for  HC  would  be  to  simulate  a  summation  of  several 
variables  in  a  control  loop.  Table  62  sl.^ws  a  flowchart  for  the  subroutine 
and  Table  188  shows  a  detailed  listing  which  includes  the  full  input/output 
list. 
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3.5.11  Subroutine  SA 

The  purpose  of  subroutine  SA  is  to  simulate  a  saturation,  A  typical 

application  would  be  to  limit  the  output  from  a  controller  to  within 

prescribed  values.  The  use  of  zero  slopes  in  the  saturation  region  should  be 

avoided  if  possible,  since  the  severe  non-linearity  will  slow  down  the  steady 
state  solver. 

Table  82  shows  a  flowchart  for  the  subroutine  and  Table  208  gives  a  detailed 
listing  which  includes  the  full  input/output  list. 

3.5.12  Subroutine  SB 

The  purpose  of  subroutine  SB  is  to  provide  the  capability  of  simulating  a 
saturation  function  with  a  dead  band. 

Table  83  shows  a  flowchart  for  the  subroutine  and  Table  209  gives  a  detailed 
listing  which  includes  the  full  input/output  list. 

3.5.13  Subroutine  Ml 

The  purpose  of  subroutine  MB  is  to  calculate  the  value  of  a  variable  according 
the  the  equation 


FO  =  C1*FNA+C2*FNB+C3*FNA*FNB 
+C4*FNA/FNB+C5 


3.5-4 


5 

I 

B 

1 


where  FO  is  the  output, 

FNA  and  FNB  are  input  variables,  and 
Cl,  C2,  C3,  C4,  and  C5  are  input  constants. 

Table  61  shows  a  flowchart  for  the  subroutine  and  Table  187  gives  a  detailed 
listing  which  includes  the  full  input/output  list. 


37 


1*- 


3.5.14  Subroutine  FV 


The  purpose  o.  ^routine  FV  is  to  provide  the  capability  of  simulating  an 
arbitrary  two-dimensional  table  Z=f{X,Y).  The  user  has  the  capability  of 
specifying  the  degree  of  interpolation  and  whether  or  not  extrapolation  should 
be  permitted. 

Table  44  shows  a  flowchart  for  the  subroutine  and  Table  170  gives  a  detailed 
listing  which  includes  the  full  input/output  list. 

3.5.15  Subroutine  RA 

The  purpose  o+  subroutine  RA  is  to  generate  random  variables  for  the  gust  wind 
model  GW.  Discrete  random  variables  with  mean  zero  and  variance  equal  to 
twice  the  time  increment  approximates  a  unit  variance  white  noise.  RA  can 
only  be  used  with  the  fixed  step  integrator  which  is  specified  by  the  command: 

INT  MODE  *  3.  A  flowchart  of  the  subroutine  is  shown  as  Table  75  and  Table 
201  gives  a  detailed  listing. 

3.5.16  Subroutine  RN 

The  purpose  of  subroutine  RN  is  to  generate  a  normally  distributed  random 
number.  This  subroutine  is  called  by  the  subroutine  RA. 

A  flowchart  of  the  subroutine  is  shown  as  Table  79  and  Table  205  gives  a 
detailed  listing. 

3.5.17  Subroutine  S2 

The  purpose  of  subroutine  S2  is  to  sum  two  sets  of  3  axis  forces  and  moments 
and  is  thus  a  summing  junction.  The  two  sets  of  forces  must  be  in  the  same 
axis  i.e.  body,  stability,  or  wind. 
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A  flowchart  of  the  subroutine  is  shown  as  Table  98  and  Table  224  gives  a 
detailed  listing. 

3.5.18  Subroutine  S3 

The  purpose  of  subroutine  S3  is  to  sum  three  sets  of  3  axis  forces  and  moments 
and  is  also  a  summing  junction.  A  flowchart  of  the  subroutine  is  shown  as 
Table  99  and  Table  225  gives  a  detailed  listing. 

3.5.19  Subroutine  TG 

The  purpose  of  subroutine  TG  is  to  transform  engine  thrust  into  body  axis 
forces  and  moments  (torques)  by  using  user  provided  direction  cosines  and 
moment  arms. 

A  flowchart  for  the  subroutine  is  shown  as  Table  110  and  Table  236  gives  a 
detailed  listing. 

3.5.20  Subroutine  TR 

The  purpose  of  subroutine  TR  is  to  transform  a  set  of  vector  quantities  from 
body  axes  to  earth  axes  using  roll,  pitch,  yaw  angle  inputs  and  direct 
multiplication  method. 

A  flowchart  for  the  subroutine  is  shown  as  Table  313  and  Table  239  gives  a 
detailed  listing. 

3.5.21  Subroutine  XP 

The  purpose  of  subroutine  XP  is  to  perform  static  transformation  on  three 
vector  quantities  (angular  rates)  from  one  coordinate  system  to  another 
coordinate  system  using  a  matrix  transformation  technique. 
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A  flowchart  for  the  subroutine  is  shown  as  Tahte  124  and  Table  250  gives  a 
detailed  listing. 

3.5.22  Subroutine  XT 

The  ourpose  of  subroutine  XT  is  to  perform  static  transformation  on  three 
vector  quantities  (torques)  from  one  coordinate  system  to  another  coordinate 
system  using  a  matrix  transformation  technique. 

A  flowchart  for  the  subroutine  is  shown  as  Table  125  and  Table  251  gives  a 
detailed  listing. 

3.6  Table  Lockup  Functions 

3.6.1  Function  TBL1 

The  purpose  of  function  TBL1  is  to  perform  linear  interpolation  on  one 
independent  variable.  It  is  a  special  purpose  function  as  it  is  called  only 
by  subroutine  TK  to  interpolate  the  trunk  shape  parameters. 

A  flowchart  for  the  function  is  shown  as  Table  102  and  Table  228  gives  a 
detailed  listing. 

3.6.2  Function  TBL2 

The  purpose  of  function  TBL2  is  to  perform  curvil inear-to-rectil inear 
transformation  and  to  linearly  interpolate  on  two  independent  variables.  This 
is  also  a  special  purpose  function  as  it  is  called  only  by  subroutines  TK  and 
AB  to  interpolate  the  trunk/bag  shape  parameters. 

A  flowchart  for  the  function  is  shown  Table  103  and  Table  229  gives  a  detailed 
listing. 
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3.6.3  Function  TBLU1 


The  purpose  of  function  TBLU1  is  to  perform  table  search  and  Lagrangian 
polynomial  interpolation  on  one  independent  variable.  The  function  uses  a 
binary  search  technique  to  locate  the  proper  position  within  a  table,  and  then 
uses  a  Lagrangian  interpolating  polynomial  of  user-defined  degree. 

A  flowchart  for  the  function  is  shown  as  Table  104  and  Table  230  gives  a 
detailed  listing. 

3.6.4  Function  TBLU2 

The  purpose  of  function  T8LU2  is  to  perform  table  search  and  Lagrangian 
polynomial  interpolation  of  user-defined  degree  on  two  independent  variables. 

A  flowchart  for  the  function  is  shown  as  Table  105  and  Table  231  gives  a 
detailed  listing. 

3.6.5  Function  TBLU3 

The  purpose  of  function  TBLU3  is  to  perform  table  search  and  Lagrangian 
polynomial  interpolation  of  user-defined  degree  on  three  independent  variables. 

A  flowchart  for  the  function  is  shown  as  Table  106  and  Table  232  gives  a 
detailed  listing. 

3.6.6  Subroutine  ETB2 

The  purpose  of  subroutine  ETB2  is  to  provide  table  look-up  for  a  function  of 
two  independent  variables  a  =  f(x,y).  The  subroutine  was  designed  to  support 
the  elastic  trunk  component  TS  and  utilizes  equal  increment  data  point  spacing 
and  linear  interpolation. 
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A  flowchart  for  the  routine  is  shown  as  Table  26  and  Table  152  gives  a 
detailed  listing. 

3.6.7  Subroutine  ETB3 

The  purpose  of  subroutine  ETB3  is  to  provide  table  look-up  for  a  function  of 
three  independent  variables  a  a  f(x,y,z).  The  subroutine  was  designed  to 
support  the  elastic  trunk  component  TS  and  utilizes  equal  increment  data  point 
spacing  and  linear  interpolation. 

A  flowchart  for  the  routine  is  shown  as  Table  27  and  Table  153  gives  a 
detailed  listing. 

3.7  Initial  Condition  Function  and  Subroutines 
3.7.1  Subroutine  IC 

The  purpose  of  subroutine  IC  is  to  solve  for  trunk  element  parametric  data  for 
free  and  loaded  ACLS  trunk  shapes  at  inital  conditions.  The  concept  is  based 
on  the  so  called  Digges  Parameter  Calculations  for  an  inelastic  trunk  (see 
AFFDl-TR -71-50,  Theory  of  an  ACLS  for  Aircraft  by  K.  H.  Digges,  June  1971). 

It  is  a  special  purpose  subroutine  called  only  by  component  TK.  The  model  can 
generate  and  store  up  to  six  sets  of  parameter  data  each  representing  an 
element  of  size  and  shape  specified  by  the  user.  Each  parameter  set  can  be 
generated  as  a  membrane  model  or  a  frozen  model  (see  Volume  I  document).  The 
trunk  section  properties  are  then  evaluated  from  the  stored  arrays  using  table 
look-up  routines  (Section  3.6)  during  the  model  execution. 

A  flowchart  of  the  subroutine  is  shown  as  Table  48  and  Table  174  gives  a 
detailed  listing. 


3.7.2  Subroutine  ICB 


The  purpose  of  subroutine  ICB  is  to  solve  for  an  air  bag  element  parametric 
data  for  free  and  loaded  air  bag  shapes  at  inital  conditions.  The  function 
solves  membrane  geometry  and  force-balance  equations  for  even  increments  of  ZO 
(the  vertical  distance  between  aircraft  hard  structure  and  bottom  of  the  air 
bag)  and  MUT  (the  bag-runway  interface  friction)  and  stores  parameter  values 
for  table  look-up  arrays.  It  is  a  special  purpose  subroutine  called  only  by 
component  AB. 

A  flowchart  of  the  subroutine  is  shown  as  Table  49  and  Table  175  gives  a 
detailed  listing. 

3.7.3  Subroutine  ICFS 

The  purpose  of  subroutine  ICFS  is  to  define  inital  condition  free  shape  Digges 
model  equations  for  an  inelastic  trunk.  It  is  a  special  purpose  subroutine 
and  called  only  by  the  iterative  equation  solver  routine  QNWT  in  subroutine  IC. 

A  flowchart  of  the  subroutine  is  shown  as  Table  50  and  Table  176  gives  a 
detailed  listing. 

3.7.4  Subroutine  ICFSB 

The  purpose  of  subroutine  ICFSB  is  to  evaluate  geometry  equations  for  a  free 
shape  air  bag  element.  The  equations  used  describe  an  inelastic  membrane 
subject  to  uniform  internal  pressure.  This  subroutine  is  called  only  by  QNWT 
in  subroutine  ICB. 

A  flowchart  of  the  subroutine  is  shown  as  Table  51  and  Table  177  gives  a 
detailed  listing. 


3.7.5  Subroutine  ICLS 


The  purpose  of  subroutine  ICLS  is  to  define  initial  condition  loaded  shape 
Digges  model  equations  for  an  inelastic  trunk.  Like  ICFS  it  is  called  only  by 
QNWT  in  subroutine  IC. 

A  flowchart  of  the  subroutine  is  shown  as  Table  52  and  Table  178  gives  a 
detailed  listing. 

3.7.6  Subroutine  ICLSB 

The  purpose  of  subroutine  ICLSB  is  to  evaluate  element  geometry  and  force 
balance  equations  for  loaded  shape  air  bag  element.  Equations  describe  an 
inelastic  membrane  subject  to  uniform  internal  pressure  and  ground 
reaction  forces.  Like  ICFSB,  this  subroutine  is  called  only  by  QNWT  in 
subroutine  1CB. 

A  flowchart  of  the  subroutine  is  shown  as  Table  53  and  Table  179  gives  a 
detailed  listing. 

3.7.7  Subroutine  KINK 

The  purpose  of  subroutine  KINK  is  to  calculate  the  arresting  gear  cable  kink 
wave  angle  at  landing  impact.  This  is  a  special  purpose  subroutine  called 
only  by  the  arresting  gear  component  AS. 

A  flowchart  of  the  subroutine  is  shown  as  Table  55  and  Table  181  gives  a 
detailed  listing. 

3.7.8  Subroutine  RES 

The  purpose  of  subroutine  RES  is  to  evaluate  the  arresting  cable  strain.  This 
is  also  a  special  purpose  subroutine  called  only  by  subroutine  KINK. 


A  flowchart  of  the  subroutine  is  shown  as  Table  77  and  Table  203  gives  a 
detailed  listing. 

3.7.5  Function  TERRA 

The  purpose  of  function  TERRA  is  to  simulate  a  rough  terrain  with  random 
profile  or  a  sinusoidal  and  (1-cosine)  bump  profiles.  For  random  profiles  the 
vertical  elevation  is  stored  in  tabular  form  as  a  function  of  position. 

A  flowchart  of  the  subroutine  is  shown  as  Table  108  and  Table  234  gives  a 
detailed  listing. 

3.7.10  Subroutine  VPRINB 

The  purpose  of  subroutine  VPRINB  is  to  print,  as  output,  the  user  supplied 
input  data  arrays  for  the  air  bag  element  dimensions,  perforations  etc.,  the 
various  stored  arrays  calculated  in  the  initial  condition  routine  ICB,  and 
element  geometry  and  forces  calculated  during  execution/simulation  as  a 
function  of  time. 

A  flowchart  of  the  subroutine  is  shown  as  Table  121  and  lable  247  gives  a 
detailed  listing. 

3.7.11  Subroutine  VPRINT 

The  purpose  of  subroutine  is  tc  print,  as  output,  the  user  supplied  input  data 
arrays  for  the  trunk  element  dimensions .perforations,  scaling  factors  etc., 
the  various  stored  arrays  calculated  in  the  inital  condition  routine  IC,  and 
element  geometry  and  forces  calculated  as  time-history  functions. 

A  flowchart  of  the  subroutine  is  shown  as  Table  122  and  Table  248  gives  a 
detailed  listing. 
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3.7.12  Subroutine  ELAS 


Subroutine  ELAS  computes  the  data  arrays  for  elastic  trunk  component  TS. 
Subroutines  ENDFS,  ENDLS,  SIDEFS  and  SIDELS  are  called  by  QNWT  to  determine 
trunk  cross-sectional  shapes.  Data  arrays  are  then  filled  with  computed  trunk 
variables. 

A  flowchart  for  the  routine  is  shown  as  Table  IS  and  Table  145  gives  a 
detailed  listing. 

3.7.13  Subroutine  ELFX 

Subroutine  ELFX  computes  incomplete  elliptic  integrals  of  the  second  kind 
E  (9,$)*  The  subroutine  was  designed  for  use  by  the  elastic  trunk  component 
(TS)  to  compute  arc  lengths  for  ellipses.  Ifftor^are  less  than  zero,  the 
routine  will  return  E  =  0. 

A  flowchart  for  the  routine  is  shown  as  Table  20  and  Table  146  gives  a 
detailed  listing. 

3.7.14  Subroutine  ELKX 

Subroutine  ELKX  computes  complete  elliptic  integrals  of  the  second  kind  E  (q, 
TT/2).  The  subroutine  was  designed  for  use  by  the  elastic  trunk  component  (TS) 
to  compute  arc  lengths  for  ellipses.  If  0  is  less  than  zero,  the  routine  will 
return  E  =  0. 

A  flowchart  for  the  routine  is  shown  as  Table  21  and  Table  147  gives  a 
detailed  listing. 

3.7.15  Subroutine  ELWR 

Subroutine  ELWR  supports  the  elastic  trunk  component  (TS)  to  provide  printout 
of  the  trunk  element  property  arrays. 
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A  flowchart  for  the  routine  is  shown  as  Table  22  and  Table  148  gives  a 
detailed  listing. 

3.7.16  Subroutine  ENDFS 

This  subroutine  is  used  by  elastic  trunk  component  TS  to  compute  the  free 
shape  of  a  trunk  end  element.  The  trunk  cross-section  is  assumed  to  be 
elliptical  in  shape  and  in  point  contact  with  the  ground.  The  shape  equations 
incorporate  membrane  theory  and  consider  elastic  deformation  in  the  trunk  hoop 
and  meridian  directions. 

A  flowchart  for  the  routine  is  shown  as  Table  23  and  Table  149  gives  a 
detailed  listing. 

3.7.1/  Subroutine  ENDLS 

This  subroutine  is  used  by  elastic  trunk  component  TS  to  compute  the  loaded 
shape  of  a  trunk  end  element.  The  trunk  cross-section  is  assumed  to  be 
comprised  of  elliptical  arcs  (2),  one  on  each  side  of  a  straight  section  which 
is  in  ground  contact.  As  in  ENDFS,  the  shape  equations  incorporate  membrane 
theory  and  consider  elastic  deformation  in  the  hoop  and  meridian  directions. 

A  flowchart  for  the  routine  is  shown  as  Table  24  and  Table  150  gives  a 
detailed  listing. 

3.7.18  Subroutine  SIDEFS 

This  subroutine  is  used  by  elastic  trunk  component  TS  to  compute  the  free 
shape  of  a  trunk  side  element.  The  trunk  cross-section  is  represented  by 
circular  arcs  and  is  assumed  to  be  in  point  contact  with  the  ground.  The 
shape  equations  incorporate  membrane  theory  and  consider  elastic  deformation 
in  the  meridian  direction  only. 
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A  flowchart  for  the  routine  is  shown  as  Table  89  and  Table  215  gives  a 
detailed  listing. 


3.7.19  Subroutine  SIDELS 

This  subroutine  is  used  by  elastic  trunk  component  TS  to  compute  the  loaded 
shape  of  a  trunk  side  element.  The  trunk  cross-section  is  represented  by  two 
circular  arcs,  one  on  each  side  of  a  straight  section  which  is  in  ground 
contact. 

A  flowchart  for  the  routine  is  shown  as  Table  90  and  Table  216  gives  a 
detailed  listing. 

3.7.20  Subroutine  XXPRT 

This  routine  is  used  by  elastic  trunk  component  TS  to  provide  printed  output 
of  trunk  variables  during  non-linear  simulation. 

A  flowchart  for  the  routine  is  shown  as  Table  126  and  Table  252  gives  a 
detailed  listing. 

3.8  Supporting  Subroutines  for  the  Foster-Hiller  Trunk  Component  (FM) 

Of  the  twenty  subroutines  which  make  up  this  section,  eighteen  were  written  by 
Foster-Miller  Associates  Inc.,  (FMA)  as  part  of  their  model  for  an  air  cushion 
landing  system  (ACLS).  The  remaining  two  subroutines  (FMWRIT  and  OUTFM)  were 
developed  by  the  Boeing  Company  to  replace  FMA  input/output  routines  which 
were  incompatible  with  the  EASY  program. 

3.8.1  Subroutine  CDVCHP 

This  subroutine  calculates  the  valve  of  DVCHP  which  is  the  ratio  of  cushion 
volume  change  to  change  in  cushion-to-trunk  pressure  ratio. 


; 
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A  flowchart  for  the  routine  is  shown  as  Table  10  and  Table  136  gives  a 
detailed  listing. 

3.3.2  Subroutine  CLRNCE 

The  purpose  of  this  subroutine  is  to  calculate  the  trunk-ground  clearance  for 
each  trunk  segment. 

A  flowchart  for  the  subroutine  is  shown  as  Table  9  and  Table  135  gives  a 
detailed  listing. 

3.8.3  Subroutine  C00RDN 

This  subroutine  calculates  the  X  and  Z  coordinates  of  the  ground  point 
corresponding  to  each  segment,  for  a  particular  ACLS  orientation. 

A  flowchart  for  the  routine  is  shown  as  Table  11  and  Table  137  gives  a 
detailed  listing. 

3.8.4  Subroutine  DYNFAN 

This  subroutine  models  the  dynamic  behavior  of  the  ACLS  fan.  Fan  pressure  is 
calculated  from  a  user  defined  polynomial  function  of  fan  flow. 

A  flowchart  for  the  routine  is  shown  as  Table  16  and  Table  142  gives  a 
detailed  listing. 

3.8.5  Subroutine  FLOW 

This  subroutine  calculates  the  ACLS  flow  and  pressure  variables  during  dynamic 
simulation. 

A  flowchart  for  the  routine  is  shown  as  Table  32  and  Table  15S  gives  a 
detailed  listing. 


49 


3.8.6  Subroutine  FMFAN 


This  subroutine  models  the  static  behavior  of  the  ACLS  fan.  Fan  flow  is 
calculated  frum  a  user  defined  polynomial  function  of  fan  pressure. 

A  flowchart  for  the  routine  is  shown  as  Table  34  and  Table  160  gives  a 
detailed  listing. 

3.8.7  Subroutine  FMWRIT 

This  subroutine  was  developed  by  the  Boeing  Company  to  print  a  list  which 
describes  each  ACLS  input  variable  and  the  corresponding  value  input  by  the 
user. 

A  flowchart  for  the  routine  is  shown  as  Table  35  and  Table  161  gives  a 
detailed  listing. 

3.8.8  Subroutine  FORCE 

This  subroutine  calculates  the  forces  and  torques  associated  with  a  particular 
ACLS  orientation. 

A  flowchart  for  the  routine  is  shown  as  Table  38  and  Table  154  gives  a 
detailed  listing. 

3.8.9  Subroutine  HYCURV 

This  subroutine  calculates  the  side  trunk  height  as  a  function  of  pressure 
ratio. 


A  flowchart  for  the  routine  is  shown  as  Table  47  and  Table  173  gives  a 
detailed  listing. 
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1.8.10  Subroutine  OUTFM 


This  subroutine  was  developed  by  the  Boeing  Company  and  replaces  FMA  routine 
OUTPUT.  The  purpose  of  the  routine  is  to  calculate  v  •.  ious  time  dependent 
variables  during  ACLS  dynami-  simulation. 

A  flowchart  for  the  routine  is  shown  as  Table  67  and  Table  193  gi /es  a 
detailed  listing. 

3.8.11  Subroutine  T ARAMS 

This  subroutine  provides  default  values  for  various  ACLS  parameters.  The  user 
may  override  these  default  values  by  inputing  the  desired  values. 

A  flowchart  for  the  routine  is  shown  as  Table  68  and  Table  194  gives  a 
detailed  listing. 

3.8.12  Subroutine  PROFILE 

This  user-supplied  subroutine  provides  values  of  ground  elevation  Yg ( i )  which 
correspond  to  ground  coordinates  Xg(i)  and  Zg(i).  A  default  subroutine  is 
included  in  the  ACLS  1 ibrarv  to  simulate  a  flat  ground  surface.  The  user  may 
input  a  different  PROFILE  using  Fortran  statements  in  the  model  description 
data  file. 

A  flowchart  for  the  routine  is  shown  as  Table  72  and  Table  198  gives  a 
detailed  listing. 

3.8.13  Subroutine  ROTATE 

This  subroutine  transforms  vector  data  from  the  vehicle  frame  to  the  inertial 
reference  frame. 

A  flowchart  for  the  routine  is  ,'h  ■  as  Table  80  ar.d  Table  206  gives  a 

detailed  listing. 
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3.8.14  Subroutine  SEGMNT 


This  subroutine  divides  the  ACLS  trunk  into  a  number  of  segments  and  assigns 
initial  properties  to  the  segments. 

A  flowchart  for  the  routine  is  shown  as  Table  84  and  Table  210  gives  a 
detailed  listing. 

3.8.15  Subroutine  SHAPE 1 

This  subroutine  provides  an  initial  assessment  of  ACLS  trunk  areas  and  volumes 
assuming  no  ground  contact. 

A  flowchart  for  the  routine  is  shown  as  Table  86  and  Table  212  gives  a 
detailed  listing. 

3.8.16  Subroutine  SHAPl 

This  subroutine  calculates  current  ACLS  trunk  areas  and  volumes  as  a  function 
of  vehicle  orientation. 

A  flowchart  for  the  routine  is  shown  as  Table  87  and  Table  213  gives  a 
detailed  listing. 

3.8.17  Subroutine  STATIC 

This  subroutine  provides  static  simulation  of  an  ACLS.  Equilibrium  conditions 
are  calculated  and  include  the  position  of  the  aircraft  center  of  gravity; 
orientation  of  the  aircraft  in  terms  of  pitch  and  roll  angles;  equilibrium 
pressure  in  the  trunk,  cushion  and  plenum;  equilibrium  flows  in  different 
parts  of  the  ACLS  and  the  theoretical  power  required  to  operate  the  ACLS. 

The  static  simulation  also  generates  load  maps  which  are  useful  in  evaluating 
behavior  of  the  ACLS  in  pure  heave  (up-down),  pitch  and  roll  modes. 


A  flowchart  for  the  routine  is  shown  as  Table  91  and  Table  217  gives  a 
detailed  listing. 

3.8.18  Subroutine  STEQU 

Subroutine  STEQU  calculates  the  time  derivatives  for  the  set  of  ACLS  state 
variables.  The  routine  calls  subroutine  FLOW  and  FORCE  to  update  the  values 
of  flows,  forces  and  toraues  associated  with  the  current  set  of  state 
variables. 

A  flowchart  for  the  routine  is  shown  as  Table  92  and  Table  218  gives  a 
detailed  listing. 

3.8.19  Subroutine  TRUNK 

This  routine  calculates  the  values  of  variables  associated  with  the  ACLS  trunk 
cross  section  shape. 

A  flowchart  for  the  routine  is  shown  as  Table  114  and  Table  240  gives  a 
detailed  listing. 

3.8.20  Subroutine  VALVE 

This  subroutine  calculates  the  pressure  relief  valve  area  as  a  function  of 
valve  displacement. 

A  flowchart  for  the  routine  is  shown  as  Table  119  and  Table  245  gives  a 
detailed  listing. 
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SECTION  IV 


COMPONENT  SUBROUTINES 


4.1  Ducting  Components 

The  component  subroutines  covered  in  this  section  include  ducts  (DU),  splits 
(FS)  and  merges  (MG). 

It  should  be  noted  that  it  is  possible  and  often  desirable  to  represent 
ducting  components  in  u*"  ar  ways.  For  example,  the  flow  merge  between  two 
bleed  air  systems  coultu  :.e  simulated  using  the  merge  component  MG  as  shown  in 
Figure  8(a).  However  if  only  one  of  the  bleed  systems  is  being  modeled  in 
detail,  and  if  the  flow  from  each  bleed  system  is  assumed  equal,  then  it  would 
be  preferable  to  use  the  multiply-and-add  standard  subroutine  MA  (see  Section 
3.5.8)  as  shown  in  Figure  8(b).  The  advantage  of  using  MA  in  this  instance  is 
that  less  input  data  is  required,  and  one  less  state  variable  would  be 
required.  Component  MA  can  of  course  also  be  used  to  simulate  a  flow  split. 
Ducts  can  often  be  omitted,  with  pressure,  temperature,  and  flow  capacitance 
effects  being  "lumped-in"  with  other  components. 

The  detailed  models  tor  ducts,  splits  and  merges  should  be  used  where  the  high 
frequency  flow  dynamics  effects  are  considered  significant;  otherwise 
alternates  should  be  used  as  described  above. 

4.1.1  Duct  DU 

The  duct  component  DU  is  a  simple  single-state  model  in  which  the  only  dynamic 
effect  simulated  is  flow  capacitance,  with  pressure  being  a  state  variable. 

The  theory  for  the  model  is  given  in  Reference  2,  Volume  I,  Section  3.1.1, 
which  also  describes  the  theory  for  more  complex  models  which  can  be  supplied 
by  the  program  user  if  desired. 
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FTgure  8  Alternate  Representations  of  Flow  Merge. 


A  flowchart  for  the  subroutine  is  shown  as  Table  14  and  Table  140  gives  a 
detailed  listing  which  includes  the  full  input/output  list.  A  warning  message 
will  be  printed  in  the  event  that  the  temperature  change  exceeds  300°F,  since 
in  this  event  a  more  complex  model  with  temperature  nodes  in  the  wall  and 
possibly  also  in  the  fluid  may  be  desirable. 

4.1.2  Split  FS 

The  split  component  FS  models  a  flow  divide,  and  has  a  single  inlet  port  and 
two  outlet  ports.  The  basic  approach  is  very  similar  to  that  of  the  duct 
component  DU.  The  flow  in  each  outlet  leg  is  calculated  based  on  input  values 
of  K  factor  and  effective  area  for  each  leg.  Heat  transfer  is  calculated 
based  on  an  input  value  of  hydraulic  diameter  (for  calculation  of  Reynolds 
number)  and  heat  transfer  area.  The  theory  for  the  Split  model  is  given  in 
Reference  2,  Volume  I,  Section  3.1.2. 

A  flowchart  for  the  subroutine  is  shown  as  Table  40  and  Table  166  gives  a 
detailed  listing  which  includes  the  full  input/output  list.  A  warning 
diagnostic  will  be  printed  in  the  event  that  the  temperature  change  exceeds 
300°F,  since  in  this  event  a  more  complex  model  with  temperature  nodes  in  the 
wall  and  possibly  also  in  the  fluid  may  be  desirable. 

4.1.3  Merge  MG 

The  merge  component  MG  models  a  flow  junction,  and  has  two  inlet  ports  and  a 
single  outlet  port.  The  basic  approach  is  very  similar  to  that  of  the  duct 
component  DU.  The  flow  in  the  outlet  leg  is  calculated  based  on  input  values 
of  K  factor  and  effective  area  for  that  leg.  Heat  transfer  is  calculated 
based  on  an  input  value  of  hydraulic  diameter  (for  calculation  of  Reynolds 
number)  and  heat  transfer  area.  The  theory  for  the  merge  model  is  given  in 
Reference  2,  Volume  I,  Section  3.1.3. 
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4'2  Fan  and  Ejector  Components 


4'2. 1  Ejector  EJ 

The  purpose  of  the  component  EJ  is  to  model  an  air  ejec-or  with  cc 
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4.2.2  Fan  with  Hysteresis  FH 


The  purpose  of  the  component  FH  is  to  model  an  externally  driven  axial  fan 
with  transition  time  constant  between  stall  and  recovery.  It  is  assumed  that 
transition  from  normal  to  stalled  operation  occurs  when  the  pressure  ratio 
exceeds  the  stall  pressure  ratio  (user  supplied)  and  that  transition  from 
stalled  to  normal  occurs  when  pressure  ratio  falls  below  a  reverse  pressure 
ratio  (user  supplied). 

A  flowchart  for  the  component  is  shown  as  Table  30  and  Table  156  gives  a 
detailed  listing. 

4.2.3  Fan  with  Surge  Analysis  FR 

The  ccmpressor/fan  model  FR  calculates  the  component  performance  using  steady 
state  maps.  Inlet  pressure  is  the  one  state  variable  for  the  model. 

The  model  is  capable  of  reverse  flow  analysis  and  can  accommodate  pressure 
ratios  less  than  1.  The  theory  for  the  model  is  shown  in  Reference  2,  Volume 
I ,  Section  3.2.2. 

A  flowchart  for  the  subroutine  is  shown  as  Table  39  and  Table  165  gives  a 
detailed  listing  which  includes  the  full  input/output  lists. 

4.2.4  ACLS  Turbofan  FT 

The  purpose  of  fan  component  FT  is  to  simulate  a  hub  or  tip  driven  axial 
turbofan  as  used  on  the  Jindivik  ACLS  vehicle.  The  subroutine  uses  input 
table  defining  steady  state  characteristics  of  tubine  and  fan.  Drive  air 
turbine  inlet  pressure  is  a  state  variable. 

A  flowchart  for  the  subroutine  is  shown  as  Table  42  and  Table  168  gives  a 
detailed  listing. 
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4.2.5  Inlet  Fan  FN 


The  purpose  of  this  component  is  to  model  a  fan  which  is  located  at  the  inlet 
of  a  duct.  It  calculates  the  fan's  flow  rate,  irlet  pressure,  inlet 
temperature  and  outlet  temperature.  The  flow  rate  is  determined  from  tabular 
data  which  is  input  as  a  function  of  pressure  ratio  (PR)  and  fan  RPM.  If  this 
flow  rate  data  has  been  normalized  for  >  where  Sis  the  air  temperature 

divided  by  the  air  temperature  at  sea  level  on  a  standard  day  and  6 is  the  air 
pressure  divided  by  the  air  pressure  at  sea  level  on  a  standard  day.  FN  will 
account  for  this  normalization  if  CORFN  is  set  to  any  value  except  0.0  and 
.99999.  The  inlet  pressure  and  temperature  are  determined  by  the  inlet  ram 
efficiency.  The  outlet  temperature  is  calcuated  from  the  fan  efficiency. 

The  user  has  the  option  of  inputing  stall  line  data  as  a  function  of  pressure 
ratio.  When  this  data  is  input,  the  program  will  print  after  each  print  cycle 
"The  fan  is  operating  in  the  stall  region"  whenever  the  fan  is  stalled. 

A  flowchart  for  the  subroutine  is  shown  in  Table  36  and  Table  162  gives  a 
detailed  listing  which  includes  the  full  input/output  lists. 

4.3  Aircraft  and  Aerodynamic  Components 
4.3.1  Generalized  6D0F  Component  SG 

The  six  degree  of  freedom  component  SG  simulates  the  aircraft  rigid  body  of 
motion  equations.  The  component  receives  as  input  all  body  axis  forces  and 
moments  including  gravity.  Parametric  inputs  include  mass,  and  moments  and 
cross  products  of  inertia  of  the  aircraft.  Outputs  from  tne  component  inciuae 
the  accelerations  and  velocities  (linear  and  angular)  in  the  body  axis 
coordinate  system,  Euler  angles  and  earth  axis  positions.  The  differential 
equations  for  the  twelve  states  representing  the  linear  and  angular  velocities 
and  positions  are  contained  in  this  component-  The  derivation  of  equations  is 
described  in  Section  2.4.1,  Volume  I.  Table  85  shows  a  flowchart  and 
Table  211  gives  a  detailed  listing. 
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4.3.2  Six  OOF  Component  DS 


This  component  is  same  as  the  component  SG  except  that  the  X-Z  plane  is 
assumed  to  be  a  plane  of  symmetry.  Therefore  the  cross  product  of  inertia 
terms  I  and  I  and  associated  terms  are  taken  out  simplifying  the  equations 

Ajr  yZ 

of  motion.  In  addition  the  airplane  positions  in  X  and  Y  earth  axis  are  not 
state  variables  and  only  rates  (velocities)  are  output.  It  is  recommended 
that  this  component  be  used  for  inflight  models  when  there  are  no  trunk  or 
ACLS  components  involved  and  the  optimal  controller  is  being  used  for  finding 
trim. 


A  flowchart  for  the  component  is  shown  as  Table  13  and  Table  139  gives  a 
detailed  listing. 

4.3.3  4DQF  Rigid  Body  Dynamics  FD 

The  four  degree  of  freedom  component  consists  of  the  same  types  of  inputs  and 
outputs  and  uses  the  same  coordinate  system  as  the  6D0F  component  DS.  The 
degrees  of  freedom  are  forward  and  side  velocities  and  roll  and  yaw  (U,V,P,R) 
and  the  model  is  applicable  to  ground  directional  control  studies.  Table  28 
shows  a  flowchart  and  Table  154  gives  a  detailed  listing, 

4.3.4  3DQF  Longitudinal  Rigid  Body  Dynamics  TL 

This  is  a  simplified  version  of  the  six  DOF  component  DS  where  only 
longitudinal  equations  of  motion  are  treated.  The  degrees  of  freedom  are 
forward  and  vertical  velocities,  and  pitch  (U,W,Q).  The  model  is  applicable 
to  long  period  longitudinal  stability  analysis  (Phugoid).  Table  112  shows  a 
flowchart  and  Table  238  gives  a  detailed  listing. 
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4.3.5  300F  Lateral  Rigid  Body  Dynamics  TO 


This  is  also  a  simplified  version  of  the  six  DOF  component  DS  where  only 
lateral  equations  of  motion  are  retained.  The  degrees  of  freedom  are  side 
velocity,  roll  and  yaw  (V,P,R).  The  model  is  applicable  to  lateral  stability 
analysis.  Table  107  shows  a  flowchart  and  Table  233  g^ves  a  detailed  listing. 

4.3.6  2D0F  Longitudinal  Rigid  Body  Dynamics  TT 

This  component  is  similar  to  the  component  TL,  above,  except  that  the  degrees 
of  freedom  are  vertical  velocity  (or  heave)  and  pitch  (W,Q).  The  model  is 
applicable  to  drop  test  simulations  and  for  inflight  short  period  longitudinal 
oscillations.  Table  116  shows  a  flowchart  for  the  subroutine  and  Table  242 
gives  a  detailed  listing. 

4.3.7  Aerodynamic  Variables  Component  VA 

The  model  adds  total  wind  perturbation  onto  the  aircraft  velocity  states, 
output  from  6D0F  component,  and  computes  aero  variables  such  as 
angle-of-attack,  sideslip,  dynamic  pressure,  and  Mach  number.  The  body  to 
stability  axis  transformation  is  achieved  wherever  applicable.  A  switch  is 
used  to  assure  compatabil ity  with  various  DOF  components.  The  derivation  of 
equations  is  described  in  Section  2.4.7,  Volume  I.  Table  118  shows  a 
flowchart  and  Table  244  gives  a  detailed  listing. 

4.3.8  Longitudinal  rorces  and  Moments  Component  0L 

Inputs  to  this  component  include  external  forces  and  torques  (e.g.  due  to 
trunk,  engine)  and  aerostabil ity  derivatives  for  the  three  longitudinal 
degrees  of  freedom.  The  external  forces  and  moments  are  combined  with  the 
linear  aero  forces  and  moments  in  body  axes  and  output  as  total  force  and 

moment  vectors  about  the  aircraft  c.g.  The  module  also  solves  for  the  linear 

• 

accelerations  U  and  W  in  order  to  compute  the  implicit  aero  terms  due  to  a  and 
8.  The  body  axis  linear  accelerations  and  the  torques  are  then  passed  to  the 
equations  of  motion  module  DS  or  SG. 
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The  aero  derivatives  can  be  specified  in  either  stability  axes  or  in  body 
axes.  For  more  details  see  Section  5.1  of  the  User's  Manual  (Reference  1)  and 
Volume  I,  Section  2.4.8.  A  flowchart  for  the  component  OL  is  shown  as  Table 
65  and  Table  191  gives  a  detailed  listing. 

4.3.9  Lateral  Forces  and  Moments  Component  DL 

This  component  is  identical  in  design  to  component  OL  and  computes  the  lateral 
set  of  forces  and  moments.  A  flowchart  for  the  component  is  shown  as  Table  12 
and  Table  138  gives  a  detailed  listing. 

4.3.10  Aerodynamic  Coefficients  Table  Lookup  AC 

This  component  provides  for  one  dimensional  table  lookup  capability  for  the 
six  basic  aerodynamic  coefficients;  lift,  drag,  pitch  moment,  roll  moment  and 
yaw  moment  as  functions  of  angle  of  attack  (a)  and  the  side  force  coefficient 
as  a  function  of  angle  of  side  slip  ($).  The  user  should  carefully  read  the 
code  before  using  this  component  or  use  FORTRAN  to  input  tabular  data  for  the 
various  non-linear  aero  coefficients.  A  flowchart  for  the  subroutine  is  shown 
as  Table  3  and  Table  129  gives  a  detailed  listing. 

4.4  Engine  and  Thruster  Components 

4.4.1  Engine  Model  (Simple)  ES 
See  4.4.2  for  description. 

4.4.2  Engine  Model  (Complex)  EC 

Two  engine  components  were  developed.  Both  calculate  engine  thrust  and  engine 
bleed-and- fan-air  pressures  and  temperatures.  The  model  ES  is  dynamically 
simple,  with  all  dynamic  effects  accounted  for  by  a  first  order  lag  for  the 
thrust  calculation. 


In  the  second,  more  complex  model,  EC,  the  engine  pressure  ratio  (EPR)  is 
computed  as  a  function  of  the  power  lever  angle.  Engine  transient  conditions 
are  accounted  for  by  introducing  a  first  order  lag  for  engine  deceleration  and 
a  second  order  lag  for  engine  acceleration.  Net  thrust  (forward  or  reverse) 
is  computed  as  a  function  of  Mach  number  and  the  actual  engine  pressure  ratio. 
The  second  order  function  requires  user  inputs  for  engine  spinup,  natural 
frequency  f(EPR)  table  and  damping  ratio  constant.  The  first  order  lag 
requires  only  the  engine  time  constant. 

The  bleed  and  fan  air  pressures  and  temperatures  are  calculated  as  functions 
of  corrected  engine  speed  which  is  a  two  dimensional  table  lookup  with  Mach 
number  and  calculated  thrust  as  independent  variables.  This  calculation  is 
common  to  both  components. 

A  flowchart  for  the  model  ES  is  shown  as  Table  25  and  Table  151  gives  the 
detailed  listing.  A  flowchart  for  the  model  EC  is  shown  as  Table  17  and 
Table  143  gives  a  detailed  listing. 

4.4.3  Yaw  Control  Thruster  YC 

The  purpose  of  component  YC  is  to  simulate  the  effect  of  an  auxiliary  yaw 
thruster  used  to  augnerrt  lateral  dynamic  stability  of  an  ACLS  equipped 
aircraft.  The  vectored  thrust  is  a  function  of  aircraft  yaw  angle  or  an  input 
control  signal.  The  vectored  thrust  may  or  may  r.ot  be  dependent  on  engine 
thrust  i.e.  the  thruster  could  be  powered  independently. 

A  flowchart  for  the  component  YC  is  shown  as  Table  127  and  Table  253  gives  a 
detailed  listing. 

4.4.4  Pitch  Control  Thruster  PT 

The  pitch  control  thruster  provides  longitudinal  dynamic  stability  to  an  ACLS 
equipped  vehicle.  The  component  is  identical  in  design  t~  YC  (4.4.3).  A 
flowchart  for  the  component  PT  is  shown  as  Table  74  and  Table  200  gives  a 
detailed  listing. 


4.4.5  Roll  Control  Thruster  RT 


The  roll  control  thruster  is  also  identical  in  design  to  the  other  two 
thrusters,  YC  and  PT,  and  provides  roll  stability.  A  flowchart  for  the 
component  is  shown  as  Table  81  and  Table  207  gives  a  detailed  listing. 

4.5  Wind  Components 

4.5.1  Gust  Wind  Model  GW 

A  wind  gus4-  model  component  developed  under  Contract  F33615-76-C-3165  has  been 
adapted  fo**  the  ACLS  library.  The  model  simulates  random  wind  gust  components 
based  on  Dryden  Spectra,  in  accordance  with  Section  3.7,  MIL-F-87858.  The 
equations  used  in  the  computer  code  are  directly  taken  from  the  MIL-F-8785B. 
Table  45  shows  a  flow  chart,  and  Table  171  gives  a  detailed  listing. 

4.5.2  Steady  or  Shear  Wind  Model  WS 

The  purpose  of  this  component  is  to  simulate  wind  shear  or  steady  wind 
components.  The  wind  magnitude  quoted  at  control  tower  altitude  of  50  feet  is 
modified  by  a  non-linear  shear  factor  to  reflect  the  change  in  wind  with 
altitude.  The  wind  vector  is  assumed  parallel  to  the  ground  plane.  The  wind 
magnitude  modified  by  shear  is  resolved  along  the  runway  coordinates  North  and 
East  (fwd  and  side)  and  transformed  into  body  axis.  The  derivation  of 
equations  is  described  in  Section  2.5.2  Volisne  I.  Table  123  shows  a  flowchart 
and  Table  249  gives  *  detailed  listing. 

4.5.3  Sunnation  of  Wind  Vectors  SV 

The  total  wind  components  are  obtained  by  the  summation  of  gust  and 
steady/shear  wind  profiles  in  the  body  axis  system  and  the  summation  is 
carried  out  in  the  component  SV.  Table  93  shows  a  flowchart  and  Table  219 
gives  a  detailed  listing. 
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4.6  Air  Cushion  System  Components 

4.6.1  Inelastic  Trunk  Model  TK 

'odel  represents  the  geometry  and  pneumatic  characteristics  of  an 
'-’j.jstic  trunk  enclosing  a  cushion  of  air.  The  trunk  is  described  by  a 
finite  number  of  elements  which  need  not  be  uniform  in  cross  section,  width, 
or  location  on  the  vehicle. 

The  trunk  is  composed  of  three  parts:  side  elements,  fore  elements  and  aft 
elements  The  side  elements  have  unconstrained  outward  movement  subject  to 
ground  contact  or  changes  in  the  cushion  „runk  pressure  ratio  and  are 
represented  with  a  membrane  model.  The  fore  and  aft  elements  experience 
little  or  no  outward  motion  due  to  constraining  peripheral  stresses  and  are 
represented  by  a  frozen  model  or  a  constrained  membrane  model  (user  option). 

The  inputs  required  are  related  to  the  trunk  geometry  such  as  section 
properties,  attach  points,  coordinates,  perforation  arrangements  etc.,  which 
can  be  obtained  -from  detailed  trunk  drawings.  The  other  data  inputs  needed 
are  trunk  damping  coefficient  as  a  function  of  trunk  flattened  area,  relief 
valve  opening  as  a  function  of  trunk  pressure  and  various  flow  discharge 
coefficients.  For  details  see  Reference  2,  Volume  I.  A  flowchart  for  the 
c  ponent  i=>  shown  as  Table  111  and  Table  237  gives  a  detailed  listing. 

4.6.2  Elastic  Trunk  Model  TS 

The  mo-'el  represents  the  geometry  and  pneumatic  characteristics  of  an  elastic 
trunk  enclosinq  a  cushion  of  air.  The  trunk  is  described  by  a  finite  number 
of  elements  which  need  not  be  uniform  in  cross  section,  width  or  location  on 
the  vehicle. 
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The  trunk  model  is  comprised  of  two  basic  units:  side  elements  and  end 
elements.  The  end  elements  take  into  consideration  membrane  loads  and 
deformations  in  both  the  hoop  and  meridian  directions.  For  the  side  elements, 
the  membrane  loads  and  strains  are  computed  for  the  meridian  direction  only; 
deformations  in  the  hoop  direction  are  assumed  to  be  zero.  As  a  user  option, 
each  element  may  be  designated  as  a  pillow  brake  element. 

Inputs  required  include  load/deflection  characteristics  for  the  membrane  as 
well  as  inputs  related  to  the  trunk  geometry  such  as  section  properties, 
membrane  attach  points,  coordinates,  perforation  arrangements,  etc.  Other 
required  data  are  trunk  damping  coefficient  as  a  function  of  trunk  flattened 
area,  relief  valve  opening  as  a  function  of  trunk  pressure,  pillow  brake 
actuation  signal,  and  various  flow  discharge  coefficients.  A  flowchart  for 
the  component  is  shown  as  Table  115  and  Table  241  gives  a  detailed  listing. 

4.6.3  Foster  Miller  Inelastic  Trunk  Model  FM  and  00 

This  model  was  developed  by  Foster  Miller  Associates  Inc.  and  adapted  by  the 
Boeing  Company  for  inclusion  into  the  EASY  program.  The  model  includes  the 
ACLS  vehicle  body  dynamics  using  two  orthogonal  coordiat->  frames  of 
reference;  an  inertial  frame  fixed  in  space  and  a  vehicle  frame  fixed  to  the 
vehicle.  The  model  also  represents  the  geometry  and  pneumatic  characteristics 
of  an  inelastic  trunk  enclosing  a  cushion  of  air. 

The  trunk  model  is  made  up  of  side  segments  and  end  segments.  Side  segments 
are  represented  by  a  simple  two-dimensional  membrane,  and  end  segments  a~e 
represented  by  a  frozen  model . 


Inputs  required  must  describe  not  only  the  trunk  geometry  as  with  component 
TK,  but  also  the  fan  performance  (static  and  dynamic)  as  polynomial 
functions,  the  ;ide  membrane  height  as  a  polynomial  function  of  pressure 
ratio,  relief  valve  dimensions,  trunk  perforation  details,  etc. 
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The  component  FM  may  be  used  in  the  Foster  Miller  mode  (user  option)  to 
duplicate  all  static  and  dynamic  analytical  capabilities  of  the  FMA/NASA  ACLS 
program.  The  component  may  also  be  used  in  the  EASY  mode  (user  option)  to 
perform  the  standard  EASY  analyses  i.e.  steady  state,  linear  analysis  and 
simulation. 

The  component  00  was  developed  to  display  variables  of  interest  during  program 
execution,  and  must  be  used  in  conjunction  with  component  FM. 

Flowcharts  for  components  FM  and  00  are  included  as  Tables  33  and  66 
respectively.  Subroutine  listings  are  provided  as  Tables  159  and  192. 

4.6.4  Air  Bag  Model  AB 

The  model  represents  the  geometry  and  pneumatic  characteristics  of  two 
parallel  inelastic  airbags.  The  bags  are  described  by  a  finite  number  of 
elements  which  need  not  be  uniform  is  cross  section,  width,  or  location  on  the 
vehicle.  The  air  bag  elements  have  unconstrained  lateral  movement  subject  to 
ground  friction  forces  in  the  Y-axis  and  are  reprcser'.ed  by  a  membrane  model. 
The  basic  design  concepts  for  the  air  bag  model  are  identical  to  those  for  the 
trunk  model.  For  derivation  of  model  equations  and  other  details,  see 
Reference  2,  Volume  I,  Section  VII. 

A  flowchart  for  the  model  is  shown  as  Table  2  and  Table  128  gives  a  detailed 
listing. 

4.6.5.  Arresting  Gear  Model  AS 

The  model  simulates  the  dynamic  response  of  a  Water  Twister  type  of  arresting 
system  composed  of  a  steel  cable  pendant,  nylon  tape  and  a  water  twister 
energy  absorber. 


rl 
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The  Water  Twister  (Registered  Trademark)  is  a  simple  water  brake  that  ^converts 
kinetic  energy  to  heat  through  turbulence.  The  brake  consists  of  a  fluid 
filled  steel  casing,  with  internal  stator  vanes,  which  houses  a  vaned 
centrifugal  rotor.  The  rotor  is  mounted  on  a  shaft  which  extends  out  of  the 
top  of  the  casing.  A  storage  reel  for  a  nylon  tape  purchase  element  is 
mounted  on  and  splined  to  the  top  end  of  the  rotor  shaft.  The  tape  is  wrapped 
on  this  storage  reel,  layer  on  layer,  forming  a  spiral  wrap.  Pulling  the  tape 
off  the  reel  causes  the  shaft  and  vaned  rotor  to  revolve  within  the  fluid 
filled  casing,  creating  turbulence.  For  derivation  of  model  equations  and 
other  details,  see  Reference  2,  Volume  I,  Section  VIII.  A  flowchart  for  the 
model  is  shown  as  Table  8  and  Table  134  gives  a  detailed  listing. 

4.7  Optimal  Controller  OC 

This  section  pertains  to  a  unique  component  in  the  EASY  library,  the  optimal 
controller  OC.  User  inputs  for  utilization  of  this  component  are  described  in 
Reference  2. 

A  complete  description  of  the  calculation  methods  and  theoretical  basis  for 
the  optimal  controller  are  presented  in  Reference  2,  Section  4.5.  Further 
elaboration  is  presented  in  Reference  2,  Volume  IV.  A  flowchart  for  the 
subroutine  is  shown  as  Table  64  and  a  listing  of  OC  is  presented  in 
Table  190. 
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4.8  Atmospheric  Flight  Conditions  FL 

The  purpose  of  model  FL  is  to  calculate  the  ambient  data  given  the  airplane's 
altitude  and  Mach  number.  Outputs  from  the  subroutine  are  the  ambient 
pressure,  ambient  temperature,  ram  pressure  and  ram  temperature.  The  ambient 
a». .  ram  temperatures  are  calculated  by  defining  the  type  of  day,  i.e.: 

Oay  =  1  MIL-STD-210B  Operation  (1  percent  risk)  Hot  Day 

2  MIL-STD-210A  Hot  Day 

3  MIL-STD-210A  Tropical  Day 

4  U.S.  Standard  Atmosphere  (Default) 

5  MIL-STD-210A  Polar  Day 

6  MIL-STD-201A  Cold  Day 

7  MIL-STD-210B  Operational  (1  percent  risk)  Cold  Day 

The  ram  pressure  is  calculated  assuming  100%  recovery. 

A  flowchart  for  this  subroutine  is  shown  in  Table  31  and  Table  157  gives  a 
detailed  listing. 
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SECTION  V 


FLOW  CHARTS 


This  section  contains  all  the  EASY  ACLS  subroutine  flow  charts  in  alphabetical 
order  by  subroutine  name. 
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Table  3:  FLOWCHART  FOR  SUBROUTINE  AC 
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Table  4:  FLOWCHART  FOR  SUBROUTINE  AF 
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Table  5:  FLOWCHART  FOR  SUBROUTINE  AMACH 
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Table  6:  FLOWCHART  FOR  SUBROUTINE  AP 


Table  6:  FLOWCHART  FOR  SUBROUTINE  AP  (CONCLUDED) 
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Table  13:  FLOWCHART  FOR  SUBROUTINE  DS 
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Table  16:  FLOWCHART  FOR  SUBROUTINE  DYNFAN 
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Table  18:  FLOWCHART  FOR  SUBROUTINE  EJ 
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Table  18:  FLOWCHART  FOR  SUBROUTINE  EJ  (CONCLUDED) 
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Table  19:  FLOWCHART  FOR  SUBROUTINE  ELAS  CONTINUED) 
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Table  19:  FLOWCHART  FOR  SUBROUTINE  ELAS  (CONTINUED) 


Table  19:  FLOWCHART  FOR  SUBROUTINE  ELAS  (CONTINUED) 
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Table  26:  FLOWCHART  FOR  SUBROUTINE  ETB2 
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Table  27:  FLOWCHART  FOR  SUBROUTIN'.  ETB3 


Table  28:  FLOWCHART  FOR  SUBROUTINE  FD  (CONCLUDED) 
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Table  30:  FLOWCHART  FOR  SUBROUTINE  FH 
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Table  30 


FLOWCHART  FOR  SUBROUTINE  FH  (CONCLUDED) 
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Table  33:  FLOWCHART  FOR  SUBROUTINE  FM  (CONTINUED) 
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Table  33:  FLOWCHART  FOR  SUBROUTINE  FM  (CONTINUED) 
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Table  33:  FLOWCHART  FOR  SUBROUTINE  FM  (CONTINUED) 


Table  34:  FLOWCHART  FOR  SUBROUTINE 
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Table  36:  FLOWCHART  FOR  SUBROUTINE  FN  (CONCLUDED) 
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Table  38:  FLOWCHART  FOR  SUBROUTINE  FORCE  (CONTINUED) 
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FLOWCHART  FOR  SUBROUTINE  FR 


Table  42:  FLOWCHART  FOR  SUBROUTINE  FT 


IF 

,IP1.  HE.cn, 


(Wl-mCAL)  /  (8S40.-FC-Y0L) 


C  CALCULATE  FAN  PERFORMANCE  FROM  INPUT  TABLES 


PSF  »  (P3-P2) -144. 

Nl-  TOT  33*4 
M2-  Tr  (23  »TOT  (3)  *4 
M3-  Tl.rra 
M4-  TOT  (3) 

M3-  60-TBLU2  (PSF.  PI . TOT  INI) .  TOT  (43 .  TOT  (M2) .  I .  I .  -M3. 
M3-AMAX1  (WtOL.HS) 

M2-  H3-H1CAL 

T3-  8«\-TlH(2»T23yV3 

T3-AMA3U  (400.. T3) 

RETURN 


-N4.N3.M4) 


169 


Table  43:  FLOWCHART  FOR  SUBROUTINE  FU 


Table  44:  FLOWCHART  FOR  SUBROUTINE  FV 
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WRITTEN  BY  -  GEORGE  DULEBfl 


INPUT/OUTPUT  LIST 


LATEST  REVISION  MAY 


TABULAR  INPUT 

_ 

INPUT  T 

OUTPUT 

ANT 

OUTPUT 

INPUT  A 

ANY 

INPUT 

TNPUT  B 

ANY 

INPUT 

DEGREE  OF  INTERPOLATION  FOR  FNA 

— 

INPUT 

A  NEGATIVE  VALUE  INDICATES  THAT  THE  NEAREST  END 
POINT  IS  TO  BE  USED  UPON  EXTRAPOLATION. 

DEGREE  OF  INTERPOLATION  FOR  FNB  -  IN 

A  NEGATIVE  VALUE  INDICATES  TiAT  THE  NEAREST  END 
POINT  IS  TO  BE  USED  UPON  EXTRAPOLATION. 


Nl=FTfl  13)  +U 
M2- FT  A  !2)  ♦FTfl  13) 

M3-FTA12) 

M4-FTAI3) 

NS=  FTA  (2)  *»ABS  !AN)  /AN 
M6-  FTA  (3)  «ABS  IBN)  /BN 
NAH*ABS  IAN) 

NBN-ABS  (BN) 

F8»TBLU2  (FNA.  FNB,  FTA  (Nl) .  FTfl  (4) ,  FTfl  (N2) .  NAN, NBN,  NS,  N6.  N3,  NU) 
RETURN 


PAGE  1 


Table  45:  FLOWCHART  FOR  SUBROUTINE  GW 


SUMUTlNC  CM  US 
l  W.KXJ.  Mt.PS.r 

x*v.x*.x#.am. 


JUS.  VS,  VSIJYS.VX.  VXS.iVX.  KS.WS. 
v^v&fslifsivsS  a».  ra.  Rxi.  ireVrs.  vs. 


Table  46:  FLOWCHART  FOR  SUBROUTINE  HI 


G™iD 

CALCULATE 
THERMAL 
CONDUCTIVITY | 
FROM  PROP 


FLOW  IN  DUCTS 
CALCULATE 
REYNOLDS  NUMBER f 
FROM  RENUX 


<  TOOOO 


LAMINAR  FLOW  OF 

AIR  INDUCTS 
CALCULATE  PRANDTL 
NUMBER  FROM  PRND 

CALCULATE 

HEAT  TRANSFER 
COEFFICIENT 

*  i 


LAMINAR 
LIQUIDS 
CALCULATE 
NUMBER  FF 

rLOW  OF 
[N  DUCTS 
PRANDTL 
iOM  PRND 

r - 

— 

- 1 

CALCULATE  BULK 
AND  WALL  VISCOSITY! 
FROM  PROP 


Table  46:  FLOWCHART  FOR  SUBROUTINE  HI  (CONTINUED) 


Table  49:  FLOWCHART  FOR  SUBROUTINE  ICS 


!W\ 

jPflGfV 


SUBROUTINE  ICS  INSET.  NPTS.A.  B.LC.GA.S1,  AMU. 
.42Q.AT3,>U.a3.AA$,|Y0.&.l.!L3,9AS) 

"T - 

JL 


PVFPOSE  -  SOLVE  FOR  AIR  SAG  ELEMENT  PARAKETERIC  DATA  FOR  FREE 

AND  LORDED  AIR  SAG  ELEMENT  SftfPES  AT  INITIAL  CONDITIONS 

METHOD  -SOLVE  ItfNIfm  GCOKETRT  AND  FORCE  BALANCE  EauATICNS 
FOR  EVEN  INCREMENTS  OF  20  AND  HUT  AND  STORES 
PARAMETER  VALUES  IN  TABLE  LOOM  IF  ARRAYS. 

LIMITATIONS  -  DUED  atY  ST  TH£  AIR  SAG  MODEL  CSFONENT  AS 
—CALL  SEOUEXCE— 

-OUTPUTS— 

ARRAYS  Y0.Ll.L3, AS  FOR  VARIOUS  20  AND  HUT 

— IBFUTS— 

NSET.NPT3.A. B.LO, GA.GB.AND  »U  tSEE  AIR  SAG  WHEN.) 


WRITTEN  BY  -  J.fi.KJLNER 


APRIL  10. 1978 


DMU-KQ/NPlS 

AMPT»*>TS 

PM-O. 

IP It) -100 


I 


H  MIC  I»l,N$gT 


PIS' 

PE33-. 

PIN) -Li  iJ 

PCn-.u.  U5-G1U) 

P 183  *.0t7U5**GA  (I) 

PCa-SIN(Pf7n 

P II 01 -CCS  IP  C7)l 

?un«siNffM8n 

p  112)  -COS  IP  (8)  1 

xii)-o.25*4.ou) 

X52)-2.0 

XC31-2.0 

JO-C 


a 


C  SOLVE  FOR  FREE  SHAPE  PAR&CTER  DATA 
C 


CAL  Q)*iT  IX.  3.5,  ICFSB.P,  .0001 .  IP,  J2.R.RHS.AJ.23 
A20 II)  -X 11) « tl.-CCS  0£  (3)  3 ) 


IF 

iRKS.LE.  .0001) 
GO  TO  IS  . 


WRITE  iitfUTE.  UU  AMS.SSCr  ! 
STOP  | 


18* 


Table  49:  FLOWCHART  FOR  SUBROUTINE  ICB 


i 


23  C9R2IUE  < 

j»»i  i 

ca  «s«  p  sta  J  •*  tie  •*  :ia 

2LI-XC3J  ; 

1XS 

«ux£u«rru-i.3-»a3*.25«sf*n.^n2J3-sw»n33)) ; 

-  'xa^xns-i.s-^as'.zs-sixa.^r.uu-jix^nici!  ; 
«0C».3«O!l!«li.-C23rii3Il)-#a^CS5»,!X53<<'i:a  i 

-  *.s«a!a-u.-3S7(;iaii**£S!  ; 

MMpAa^  is  •*»  as  *.  w  a  *  rs 


IF 
n.Ea.a 


(c  nag  nsnuf.  m  i»  a  rnrnn 


L 


mu.  K.n«rra! 

- i*«.i ! 


«u.!.S*W; 


*<<  i 


. . . ««»««•»» 


fSSS 


Table  49:  FLOWCHART  FOR  SUBROUTINE  ICB  (CONTINUED) 


(c  store  negative  m  m  i  araat s') 


nou.K.n*mro 

K.1  U,K.I1*X1 
IL3U,K,I)»X3 
*A3U,K.I)«DA3 


s: 


CONTINUE 


(iHAtTC.eooa  i 

uhaite.  soon  nnts.wu.a  in  ,i  in  j.o  in .  ga  id  ,  g*  in 

tlMAITE . 800S1 AM (1)  .ATOU.  1.11  ,AU  II.  l.il  .AftSU.  1.11 
ilMUTE.OOtO 


X 


L{ 


H  88  10  K»l. NATS  I 

TI 


«  o;  oo  *«o  in-'urrwcMi  /anpt  I 
- 1 - 


WnEUWIITE,  60201  (M  00  .K»l. IPT3> 
WAIT?  (IMAIT?,  8030 


X 


- - M  SO  >t.NP73| 

I  Nwmiii»inE'.602ffl  wTou.Kln.Kxt.iyTa  ) 


-fso 


CONTINUE 

I 


I  NAflE  IIHRITE.  SOS) 

'-JSSLX - 


BO  SI  >1.  NATS  | 


_ _ X _ , 

I  wuTEawuTE.acam  iirou.k.n.oi.irra  ] 

- (El  CONTINUE  1 


PAGE 

SCJ 
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Table  49:  FLOWCHART  FOR  SUBROUTINE  ICO  (CONCLUDED) 


rg 


r* - »  si  >t.*rT3  j 

j  — - 1 

j  *K\  u.r,.  n  .sTi. >r^i  I 

5  X 


p*fu  awtTt.eoauT! 


, - , - -fw  si  j-* ,  fcr  tjj 

I  1.  *nEiivnTfTi?ro3)  v*,3 u.tCliTKM.y*?!  | 

\ _  - — I L 


cwriwun 


»wm  iiwrTiwFI 

- r— - 


-H  »  a  >r^n 
i 


5  [  WTttiNWnt  60W1  UMU.K.  1)  .KMTyrTn 

L  - L 


HniTtuwm.cewj 

j=fczr 

~W  N?  SO  >\,W»T3  I 

X 


L. 


1  ^ITtUHUtlC.vwa  g^su;^.  !),!<» \,yTS~! 


I 


— (m  cflStwTl 

,_z£_ 

1  mtTtuwnt.tiwa  1 

~W  S',  J-i.frT;  [ 


I  (..wttuwnt.wga  (W}u,n, n ,k»uh*ts) 


r  & 


Table  50:  FLOWCHART  FOR  SUBROUTINE  ICFS 


SUBROUTINE  ICFS(X,M,K,R.P) 


PURPOSE  -  DEFINE  INITIAL  CONDITION  FREE  SHAPE  DIGGE5  MODEL 
-  FAR  ACLS  TRUNKS 


MEHTOB  -  SEE  AFFDL-TR-71-S0.  THEORY  OF  AN  ACLS  FOR  AIRCRAFT 
-  BY  KENNERI.Y  H.  DIGGES.  JUNE  1971. 


LIMITATIONS  -  USED  ONLY  BY  ROUTINE  QNHT  IN  SUBROUTINE  IC 


WRITTEN  BY  -  J.R.KILNEH  AND  H.K.WAMI 


OCT.  21  1977 


Am-xm-xc2i«a.-P(5n 

R  (23  *P  'U)  -X 11)  *<X  O)  -X  (23  "X  lU) 

a  raj -p  (2) -x  in  «sin  ix  an -x  (2) -sin  a  wn 

R  N)  (3)  -X 11)  -  U .  -CCS  (X  (3)) )  «*  (2) « (1 .  -COS  (X  cun ) 
RETURN 
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ICFS 
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Table  51:  FLOWCHART  FOR  SUBROUTINE  ICFSB 


\  SUBROUTINE  ICFSB  tX.N.K.R.PJ  / 
- - 


PURPOSE  -  EVALUATE  GEOMETRY  EQUATIONS  FOR  FREE  SHAPE  AIR  BAG 

-  ELEMENT 

METHOD  -  EQUATIONS  DESCRIBE  INELASTIC  MEMBRANE  SUBJECT  TO 

-  UNIFORM  INTERNAL  PRESSURE 

LIMITATIONS  -  USED  ONLY  BY  ROUTINE  QNWT  IN  SUBROUTINE  ICB 
u  h  tt  OUTPUT  h  m 

X  THE  SOLUTION  ESTIMATE  FROM  THE  LAST  ITERATION 

R  THE  RESIDUAL  VECTOR 

<MM<INPUTS<MW 

X  INITIAL  ESTIMATE  VECTOR 

Pit)  PRINT  CONTROL  OPTION 

P-0  NO  PRINT  OUT 

P-K  FOR  PRINT  OUT  AT  EVERT  K-TN  ITERATION 
P 12) . .  ,P  (N)  STORAGE  AVAILABLE  FOR  PASSING  DATA 
TO  SUBROUTINE  FUNIICFS)  BY  QNWT. 

M.K  DIMENSION  INDICATORS  FOR  X  AND  AJ  MATRICES 


WRITTEN  IT  J.R.KILNER 


APRIL  10.1978 


R(l)*PM)-X(l)«(X(2)«-X(3)) 
fl  (2)  -P 12)  -X  m  « (SIN  002) )  <-SI.V  003) ) ) 
fl  13)  f3)  -X  (1) « (COS  003) )  -COS  tX  f2) )  5 

RETURN 


ISO 
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Table  52:  FLOWCHART  FOR  SUBROUTINE  ICLS 


\  SUBROUTINE  ICLStX.H.K.fi.PI  / 


PURPOSE  -  DEFINE  INITIAL  CONDITION  LORDED  SHAPES  DIGGES  MODEL 

-  FOR  RCLS  TRUNKS 

MEHTOD  -  SEE  HFFDL-TR-71-SO.  THEORY1  OF  RN  RCLS  FOR  RIRCfifiFT 

-  ST  KENNERLY  H.  DIGGES.  JUNE  1371. 

LIMITATIONS  -  USED  ONLY*  BY"  ROUTINE  QNHT  IN  SUBROUTINE  IC 
WRITTEN  BY  -  J.R.KILNER  AND  M.K.NRHI  OCT.  21  1377 

_ y 


RU)-XU)-X<2JwU.-*»Cn 

R  (2J  PU  -X  U)  -X  C3J  -X 12)  *X  N)  -X  CS3 

Roi  -e  t2J  -yu)  msim  ccoi )  -xgj  -sin  or  «)  i  -xra 

R  Ml  «P  £31  -X lit  *r  11 . -COS  003) ) )  *X O  •»  11  .■ -COS  Cf  M) ) ) 

R  IS -P 1EJ -X  OI  f  ll.-COS  Of  M)  n 

RETURN 

191 
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Table  53:  FLOWCHART  FOR  SUBROUTINE  ICLSB 


\  SUBROUTINE  IO.SB1X.N.K.A.P)  / 

_  '""T ^HL _ 

PURPOSE  -  EVALUATE  ELEMENT  GEOMETRY  AND  FORCE  BALANCE  EQUATIONS 
-  FOR  LORDED  SHAPE  AIR  SAG  ELEMENT 


HEWm 

LIMITATIONS  -  USED  ONLY  BY  ROUTINE  QNHT  IN  SUBROUTINE  ICB 

FOR  INPUTS  AND  OUTPUTS.  SEE  o’iOOUTINE  ICFS 
WRITTEN  BY  -  J.R.KILNEfi  APRIL  10.1978 


- >4  DO  10  1-1.2 

I  ZE_ _ 

p  uniu  -x  u*2i  -x  in  - is.  mt6HMi»6)) 


^  IF 
•JP  UH*l)  .LT.O.L 


pih-i2)=*3.iui6*pii»6) 


ptmui-o. 


10  IF 
•UnU)  .LE.O.L 


PII*12)  »X  (1*2) /XII) 


Rti)«xti)-xt2)-P(S)-xis} 

R  (2)  -P  (4)  -X  (3)  -X  ?4)  -X 151 
R  (3)  -P  (2)  -X 115  -SIN  IP  113) )  -X 12)  -SIN  IP  11U) )  -X  IS) 

-  *P  US)  *#*  110)  *P 116)  -P  (12) 

R  (U)  «P  (3)  -X 11)  - 11 .  -COS  (P 113) ) )  ♦X  (2)  - 11. -COS  tP  UU) ) ) 

-  ♦PUS) (9) -P 118) -PUD  „„„ 

R  (S)  -P  (8)  -X  (2)  -  (l .-COS  IP  UR) ) )  ♦?  11 8)  -P  ll  1) 

RETURN  _ _ 
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Table  54:  FLOWCHART  FOR  SUBROUTINE  IT 


\  SUBROUTINE  inFO.FODOT.IFO.FIN.GKI.GKL.RMA.AMI)  / 

. . .  I 

VERSION  2.  REVISED*  OCT  8  1976 

PURPOSE  -  simulation  OF  an  integrator  with  saturation 

METHOD  -  SEE  CODING 

LIMITATIONS  -  EXCESSIVELY  HIGH  VALUES  OF  GKL  MAY  RESULT  IN  POOR 
STEADY  STATE  CONVERGENCE 


HRITTEN 

BY  -  ADAM  LLOYD 

LATEST  REVISION  -  NOV  IS 

INPUT/OUTPUT  LIST 

FO 

INTEGRATOR  OUTPUT 

ANY 

OUTPUT  STATE 

FQD&T 

OUTPUT  DERIVATIVE 

ANY 

OUTPUT  DERIV 

IFO 

INTEGRATOR  CONTROL 

— 

PROGRAM  VAR 

FIN 

FUNCTION  INPUT 

ANY 

INPUT  VAR 

GKI 

INTEGRATOR  GAIN 

ANY 

INPUT  PARAN 

GKL 

DERIVATIVE  LIMITER  GAIN 

ANY 

INPUT  PRRAM 

AMA 

UPPER  LIMIT  OF  OUTPUT 

ANY 

INPUT  PFRAM 

WHERE  DERIV.  LIMITER  STARTS 

AMI 

LOWER  LIMIT  OF  OUTPUT 

ANY 

INPUT  PARAM 

CL 


I 


WERE  DEAIV.  LIMITER  STARTs') 
- 


EPS=FIN 


SPACE  2 
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Table  55:  FLOWCHART  FOR  SUBROUTINE  KINK 


\  SUBROUTINE  KIWtVO.CSi.EC.DNC.PKH.C.X)'/7 


PURPOSE  -  DETERMINE  ARRESTING  CABLE  KINK  HAVE  ANGLE 

MEHTOD  -  SEE  ARRESTING  SYSTEM  DOCUMENTATION. VOL  I.  FINAL  REPORT 

LIMITATIONS  -  CALLED  OftJ  BT  THE  ARRESTING  SYSTEM  COMPONENT  AS 


CALL  SEQUENCE 
INPUTS  mm. 


VO 

CSI 

EC 

DNC 


VEHICLE  VELOCITY  .  FT/SEC 
INITIAL  CABLE  STRESS,  LBS/SG  IN 
CABLE  MODULUS  OF  ELASTICITY,  LSS/SQ  IN 
CABLE  HEIGHT  DENSITY.  LBS/CU  IN 


OUTPUT 


PKH  -  KINK  HAVE  ANGLE.  RAD 

C  -  CAM.E  SONIC  VaCCITY.  FT/SEC 

CSTR  -  STRAIN  IN  HAKE  OF  TRANSVERSE  HAVE.  IK/IN 


WITTEN  BT 


J.fl.KILNER 


MAT  1378 


OSQflT  (3M.«EC/DNO/l2. 

VOVO/C 

Pt2J-VC«VC 

PCSJ-CSI/EC 

XU5-Pt3) 

P  (U  *0. 

IPUJ*SO 

CRLL  QNHTOC.l.l.RES.P.  .0001,  IP. 0,  A. RMS. AJ.  SB 
SINPKH-VC/SWT  a  Cl)  -  tl .  *X  U) » 
PKH-ASINtSINPKW 
RETURN 


195 
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Table  56:  FLOWCHART  FOR  SUBROUTINE  LA 


\  SUBROUTINE  LAlF8.FaCflT.IF0.FIM.GflI. TO  / 

PURPOSE  -  TO  SIMULATE  FIRST  ORDER  LAG  FO  GAI 

FIN  11 .  *TC«*SY 

METHOD  -  SEE  CODING 

WRITTEN  5T  -  ADAM  LLOYD  LATEST  AEVISJON  NOV  7S 

LIMITATIONS  -  TIME  CONSTANT  TC  SHOULD  NOT  EQUAL  ZERO 


INPUT/OUTPUT  LIST 


FO 

TRANSFER  FUNCTION  OUTPUT 

ANY' 

OUTPUT  STATE 

FODCT 

TRANSFER  FUNCTION  OUTPUT  DERIV. 

ANT 

OUTPUT  STATE 

IFO 

1MTEGERAT0R  CONTROL 

— 

PROGRAM  VAR 

FIN 

TRANSFER  FUNCTION  IFPUT 

ANT 

INPUT  VAR 

GAI 

TRANSFER  FUNCTION  GAIN 

— — 

INPUT  PARAH 

TC 

TIME  CONSTANT 

SECS 

INPUT  .’ARAN 

196 
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Table  59:  FLOWCHART  FOf.  SUBROUTINE  LL 


C 

c 

r 

C 

c 

c 

c 

c 

r 

C 

C 

r- 

U 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ggi.xim  IXl.FO.FIW.TCl.TC2.Gair7 

“PURPOSE  -  TO  SIMUTTE  LEAD  LAC  TRANSFER  FUNCTION - 

FC  GAI«  tl.*TCl«S) 

FIN  ll.-TC2«$) 

METHOD  -  SELF  EXPLANATORY 


LIMITATIONS 


NONE 


witten  sr  -  a  nan  lur’d 

input/output  list 
xi 

XI  DOT 
IXl 
FO 
FIN 
TCI 
T C2 
GAI 


LATEST  REVISION  NOV  75 


STATE  VARIABLE 
STNTE'SWflJLE  DERIVATIVE 
INTEGRATOR  CONTROL 
TRANSFER  FUNCTION  OUTPUT 
TRANSFER  FUNCTION  imJT 
TIME  CONSTANT  [NUMERATOR) 
TIME  CONSTANT  [DENOMINATOR) 
TRANSFER  FUNCTION  GAIH 


m 

ANT 

ANT 

ANT 

SECS 

SECS 


T 


OUTPUT  STATE 
OUTPUT  STATE 
PROGRAM  VAR 
OUTPUT  VAR 
INPUT  VAR 
INPUT  PfiRAM 
INPUT  PARAM 
INPUT  PARAM 


FO=  DCI»FIN«TC1«GAI)  /idT] 


f  RETURN  j 


199 


PAGE  1 


no noon non nnnnnnnnn 


Taole  50: 


"LOW CHART 


FOR  SUBROUTINE  MA 


V  SJ3R3UTINE  nfi  IF3.FI".C1.  C22  x 

-  — , - 

J. 


y 


PURPOSE  -  13  SIMULATE  THE  ESUfiTiaN  0UT?UT«Cl«ifPUT 


METHOD  -  SEE  CODING 

BITTEN  3T  -  ADAM  LL3T3  LATEST  REVISION 

LIMITATIONS  -  NONE 


INPUTyflUTPUT  LIST 


F0 

OUTPUT  VARIA5LE 

FIN 

INPUT  VARIABLE 

F&X 

Cl 

CONSTANT  MULTIPLIER 

— 

C2 

CONSTANT  ADDITION 

— 

!  FC*Cl«f iN  -  C2 
1  FgTURN 


200 


OUTPUT  VAR 
INPUT  yar 
INPUT  PflRflK 
iSPUT  PARfSl 
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Table  62:  FLOWCHART  FOR  SUBROUTINE  MC 


SUBROUTINE  HC  iFa,FlN,F10.FIP,C\,C2,C3.Cm  / 


X~E - PUfVOSE  -  TO  SIMULR'tTThE  EQUATION  F0*Ct**FlN*-C2«F10»C3«nP*'CH 

/  C 

c  METHOD  -  SEE  CODING 


WHITTEN  BY 


ADM  LLOYD 


LATEST  REVISION  NOV  75 


LIMITATIONS  -  NONE 


INPUT/OUTPUT  LIST 


OUTPUT  VARIABLE 
INPUT  VARIABLE 
INPUT  VARIABLE 
INPUT  VARIABLE 
CONSTANT  NULTIPUER 
CONSTANT  MULTIPLIER 
CONSTANT  MULTIPLIER 
CONSTANT  ADDITION 


OUTPUT  VAR 
INPUT  VAR 
INPUT  VAR 
INPUT  VAR 
INPUT  param 
INPUT  PARAM 
input  PARAM 
INPUT  PARAM 


F0*'Cl«FlN»C2wn0*C3*<nP»CU  j 
RETURN  _ ! 
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Table  63:  FLOWCHART  FOR  SUBROUTINE  MG 


?03 


Table  64:  FLOWCHART  FOR  SUBROUTINE  OC 


\  SUWOUTINE  OClG.S.AK.fN.X.XDOT.  INT.r.rOP.U,UOP.N;Nr.NU]  / 

zzzzx^ 


VERSION  l. 


Rfrflsc t  optiiwl  controller  component 

CALL  SEflUENCEi  0  -  OUTPUT  MATRIX 


REVISED)  MAT  0  197S 


s  -  input  matrix 

*  -  STRIILITT  MATRIX 

n<  -  B.C.  SAIN  MATRIX 

X  -  STATE  VECTOR 

XDOT  -  RATE  VECTOR 


iMT 

r 

TOP 

u<jp 

N 
NT 

NU  _ 

cesiOMED  rrt  j. d. surrough? 

—  CALCULATE  O.C.  RATES 


-  INTEGRATOR  CONTROL  VECTOR 

-  INPUT  VECTOR 

-  INPUT  OPERATING  POINT  OPPSET  VECTOR 

-  OUTPUT  VECTOR 

-  OUTPUT  OrtJVrTING  POINT  OfESET  VECTOR 

-  OPTIMAL  CONTROLLER  ORDER 


NLf<|£A  Of  INPUT^ 


-  NUMBER  Of 


TEST  EOR  ZERO  ORDER  CONTROLLER 


APRIL  1373 
XDOT  •  S«T*-AK«X 


CALCULATE  OUTPUTS 

‘ '  3 


H»o  »  iw  i«iTnu) 

|  SJMaJOP  in  1 

X 


??o  >i,nt  1 


X 
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Table  64:  FLOWCHART  FOR  SU3R0UTINE  OC  (CONCLUDED) 


IfWE  1  | 
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Table  65:  FLOWCHART  FOR  SUBROUTINE  OL 


Table  67:  FLOWCHART  FOR  SUBROUTINE  OUTFM 


ioa  amnu 


Table  67:  FLOWCHART  FOR 


SUBROUTINE  OUTFM  (CONTINUED) 


Table  67:  FLOWCHART  FOR  SUBROUTINE  OUTFM  (CONCLUDED) 


lib  CTim'l 


r?yFit>TniTrrf;fii!rmr 


Table  74:  FLOWCHART  FOR  SUBROUTINE  PT  (CONCLUDED) 


PflCt  ! 

PT 


219 
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Table  75:  FLOWCHART  FOR  SUBROUTINE  RA 


\  SUBROUTINE  RA  tXNU.XNV.XNN.XNP)  /' 


VERSION  2.  REVISED  OCT  22  1976 

PURPOSE*  GENERATE  RANDOM  VARIABLES  FOR  WIND  MODEL 
METHOD*  DISCRETE  RANDOM  VARIABLES  WITH  MEAN  2ER0  AND 

VARIANCE  »2«TINC  APPROXIMATE  UNIT  VARIANCE  WHITE  NOISE 


OUTPUTS 

XNU.XNV.XNN  -NOISE  SAMPLES  FOR  U.V.W  GUST  VELOCITIES 
XNP  -NOISE  SAMPLE  FOR  P  ANGULAR  RATE  GUST 

DESIGNED  BT  A. H. WARREN  SEPT.  1275 


SIG-  SORT  mNC^TiNCl 
CALL  RN IXNU.  DUN.  SIG.  0.3 
CALL  RNOCNV.DUH.SIG.O.) 
CALL  RNKNW.DUM,  SIG. 0.3 
CALL  RN  (XNP,  DUM,  SIG.O.) 
RETURN 
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Table  77:  FLOWCHART  FOR  SUBROUTINE  RES 
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Table  79:  FLOWCHART  FOR  SUBROUTINE  RN 


\  SUBROUTINE  ANOJ.AX.SIG.AMW  / 
- 1 - 


VERSION  l.  REVISED.  OCT  7  1976 

PURPOSE  -  GENERATES  A  NORfWxf  DISTRIBUTED  RANDOM  NUMBER 

ra|  i  ccn*|c*i<-£r 

U  -  the  RANDOM  NUMBER  OUTPUT 

AX  -  A  START  PARAMETER  WHICH  CONTROLS  THE  BEGINNING  POINT 
OF  THE  OUTPUT  SEQUENCE.  AX  SHOULD  RE  ANT  ODD  NUMBER 
GREATER  THAN  ONE.  THE  DEFAULT  VALUE  Or  AX  IS  N3W6971. 
AX  IS  UPDATED  FOR  NEW  CALLS  TO  THE  SUBROUTINE 
SIG-  THE  DESIRED  STANDARD  DEVIATION  OF  THE  SEQUENCE 
AMN-  THE  DESIRED  MEAN  OF  THE  SEQUENCE 
DESIGNED  BY  ROGER  H.  CALL  SEPT.  1976 


IX-IXkIY 

SU1*SUM*‘IX/281U7U97671O0SS. 


-|l  AX-IX 


U-tSUM-6.ffl«SIG*AMN 

RETURN 
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5WSJWC: 


Table  82:  FLOWCHART  FOR  SUBROUTINE  SA 


SUEROUTINE  SA(FQ,FIN,C1,C2,C3,CU,CS,C6) 


PURPOSE  -  TO  SIMULATE  SATURATION 


METHOD  -  SEE  CODING.  C3  AND  C8  ARE  VALUES  OF  THE  INPUT  AT  WHICH 
SATURATION  OCCURS.  C3  IS  GREATER  THAN  C8.  THE  ROUTINE 
CAN  SIMULATE  A  CHANGE  OF  SLOPE  AT  THE  ORIGIN  (C1.NE.C4) 
PROVIDED  C6  IS  LESS  THAN  ZERO.  SIMILARU  THE  SLOPES 
IN  THE  SATURATION  REGION  (C2  AND  CS)  CAN  DIFFER. 

THE  SLOPES  CAN  BE  POSITIVE  OR  NEGATIVE 


WRITTEN  BY  -  ADAM  LLOYD 


LATEST  REVISION  -  NOV  75 


LIMITATIONS  -  USE  OF  ZERO  SLOPES  IC2-0  OR  CS-O)  IN  THE  SATURATION 
REGION  SHOULD  BE  AVOIDED.  IT  IS  DESIRABLE  THAT  THE 
SLOPE  RATIOS  C1/C2  AND  C4/CS  SHOULD  NOT  EXCEED  100. 
EXCESSIVE  SLOPE  RATIOS  MAY  RESULT  IN  VERY  5L0W 
CONVERGENCE 


INPUT/OUTPUT  LIST 


OUlrUT  VARIABLE 
INPUT  VARIABLE 
SLOPE 


)  FIRST 


OUTPUT  VAR 
INPUT  VAR 
INPUT  PARAM 


_ ] _ 

^c 

C3 

SATURATION  INTERCEPT) 

ANY 

INPUT 

PARAM  ^ 

c 

CU 

SLOPE  ) 

SECOND 

ANY 

INPUT 

PARAM 

c 

CS 

SAITJRATION  SLOPE  ) 

SLOPE 

ANY 

INPUT 

PARAM 

c 

CP 

SATURATION  INTERCEPT) 

ANY 

INPUT 

PARAM 

IF  ^ 
'FIN  oT.C3) 

^  GO  TO  10 


IF 

(FIti.LT.C0) 
^GO  TO  20  . 


if 

(FIN.LT.O.) 
^  GO  TO  30  ^ 


C  NEGATIVE  UNSATURATED 


FO»Cl*FIN 
GO  TO  100 


(C" 'nTgATIVE  SATURATION 


30  FO-CUkFJN 


100  RETURN 


20  F0-CU«C6'CS»«(FIN-C6) 

GO  TO  100 
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Table  83:  FLOWCHART  FOR  SUBROUTINE  SB 


\  subroutine  sr(fo.nN.cT^C2.c3.cw.csTc8~.c7.cs)  y1 


WttE  -  TO  SIMUJITE  SATURATION  WITH  A  OCAO  3AN5 


KETH6B  -  SEE  COOIMG. 


tHITTEK  IT  *  GCCRGC  3ULEM 

LIMITATIONS  -  Cl  MUST  If  POSITIVE 

C3  *VST  AC  GREATER  IRAN  Cl 

cs  *}sr  At  positive 

C7  fiutl  K  GREATER  THAN  0. 
r.  MUST  BC  NEGATIVE 

d  must  sc  less  turn  C2 

CS  MUST  BE  NEGATIVE 
C3  MUST  AC  GREATER  THAN  0. 


latest  rcvisiom  -  fes  it 


iwvi/ainPUT  list 

FC 

SJV'ft  VflRIHBU 

ANT 

OUTPUT  VAR 

FIH 

IJ*U?  VARIABLE 

aht 

jNPUT 

VAR 

Cl 

POSITIVE  DEAD  BAND  ON  FIH 

ANT 

l.Stf*UT 

VAR 

C2 

NEGATIVE  DEAD  BANJ  OH  FIN 

ANT 

Itf’VT 

VhR 

a 

POSITIVE  SATURATION  INTERCEPT 

ANT 

INPUT 

VAR 

CH 

NEGATIVE  SATURATION  INTERCEPT 

ANT 

IHVT 

VAR 

CS 

POSITIVE  SATURAT ION  LIMIT  OH  FO 

ANT 

INPUT 

VAR 

CS 

NEGATIVE  SATURATION  LIMIT  CM  FO 

ANT 

INPUT 

VAR 

C7 

SATURATION  SLOPE 

ANT 

INPUT 

VAR 

C3 

SATURATION  SLOPE 

ANT 

INPUT 

VAR 

TP»  SUER0-C1 
TN»  SUERO-tt 
SU€G-  (C0-TN1/  CCU-C2J 
SLPLUS-  (CS-TPJ  /  (C3-C1J 
FO  «SL2EP0-fIW 


Table  84:  FLOWCHART  FOR  SUBROUTINE  SEGMNT  (CONTINUED) 


Table  84:  FLOWCHART  FOR  SUBROUTINE  SEGMNT  (CONCLUDED) 


JWFS'CfFIV 


FLOWCHART  FOR  SUBROUTINE  SHAPE  2  (CONTINUED) 


Table  87: 


(cT5~ carnet  srwuaff  xcTipaT 


flowchart  for  subroutine  shape  2  (continued) 


Table  87:  FLOWCHART  FOR  SUBROUTINE  SHAPE  2  (CONTINUED) 


Table  87:  FLOWCHART  FOR  SUBROUTINE  SHAPE  2  (CONTINUED) 
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Table  87:  FLOWCHART  FOR  SUBROUTINE  SHAPE  2  (CONTINUED) 


W3 


can  mm  l 


(c  camitt  rium  m  omx) 
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Table  87:  FLOWCHART  FOR  SUBROUTINE  SHAPE  2  (CONTINUED) 
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Table  91 : 


FLOWCHART  FOR  SUBROUTINE  ST AT I L 
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Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 
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Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


mmssL 
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Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 
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Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


Table  91 


:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 
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TOTIHB  VI 
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Table  91:  FLOWCHART  FOR  SBURGUTINE  STATIC  (CONTINUED) 


Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


r 


91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 
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Table  91:  FLOWCHART  FOP  SUBROUTINE  STATIC  (CONTINUED) 


net  ii 


Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


na  n 


Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


me  h 


Table  91 : 


FLOWCHART  FOR  SUBROUTINE  STATIC  (CONTINUED) 


STATIC  (CONTINUED 


f  oa  tal 


ms  » 

STATIC 
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Table  91:  FLOWCHART  FOR  SUBROUTINE  STATIC  (CONCLUDED) 


X 


pm  * 

STATIC 
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Table  92:  FLOWCHART  FOR  SUBROUTINE  STEQU 


\  xmanm  sttw/ 


stmt  nutria*  rm  m  rm wc  ststoi 


fttLwiNS  imm«  m  ollo  to 
Ftfca.Tmajo  «•  aaa.  stva  tx  mu 
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9  TX 
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\*srrm71m  aw  itemwz  fum 

■mmiwa  in  tie  cwnao  m  urn  to  accatamt  artci*.  caamws 


ww.  sn 


(c  catnap  Ttutt-wjgm  grwwia  .  utmiI 


to  wtrm-o»<«mj*aoTi/najtxTio««r»eMPUJwnjn'-«TO*-«Twrrl 
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FOR  SUBROUTINE  STEQU  (CONCLUDED) 


Table  93:  FLOWCHART  FOR  SUBROUTINE  SV 


SUBROUTINE  $V  OJH,  VW,  WH,  PH,  OH,  RH.  UHS,  VHS,  HHS,  UG,  VG,  WQ,  PG,  QG, RG)  7 


X 


IS 


l.  NAY  26  1977 

SUN  TWO  SETS  OF  3  AXIS  VELOCITIES  AND  ANGULAR  RATES 
METHOD  ADD  STEADY  OR  SHEAR  WIND  COMPONENTS  TO  THE  RANDOM 
GUST  COMPONENTS  AND  GUST  ANGULAR  RATES 
CALL  SEQUENCE 
OUTPUT  3 

LINEAR  VELOCITIES  —  BODY  AXES 
UN.VW.WM  -SUM  OF  X.Y.Z  AXIS  HIND  VELOCITIES,  FT/SEC 
ANGULAR  VELOCITIES  —  BODY  AXES 
PH.QH.RH  -SUM  Or  X.Y.Z  AXIS  ANGULAR  VELOCITIES,  DEG/SEC 
hkwwk  I HPUT  S 

.LINEAR  V&OCITIES  —  BODY  AXES 

UHS.VHS.HHS  -X.Y.Z  AXIS  STtADT/SHEAR  HIND  CSMPONENTS,  FT/SEC 
UG.VG.HG  -X.Y.Z  AXIS  GUST  HIND  COMPONENTS,  FT/SEC 
ANGULAR  VELOCITIES  —  BODY  AXES 
PG.QG.RG  -X.Y.Z  AXIS  GUST  ANGULAR  COMPONENTS,  DEG/SEC 
MRImEN  BY  MAHINDOT  H&I  MAY  1977 

— «  SUM  LINEAR  VELOCITIES  «*♦* 


UW»  UM$  -  UG 
VW-  VN3  «•  VG 
HW-  WHS  WG 


d 


SUN  ANGULAR  VELOCITIES 


D 
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Table  94:  FLOWCHART  FOR  SUBROUTINE  SW 


SW  -  ABS 


Table  95:  FLOWCHART  FOR  SUBROUTINE  SX 


SUBROUTINE  SX  IVOl , V02. VAl , VA2, VB1, VB2, SHI . TCI , TC2) 


X  c 
/  c 

PURPOSE  - 

TO  PROVIDE  A  SWITCH  COMPONENT 

FOR  TWO  VARIABLES 

/  C 

/  C 

/  c 

METHOD  - 

SEE  CODING 

/  c 

/  c 

/  c 

WRITTEN  BY 

-  ADAH  LLOYD  LATEST  REVISION 

NOV  7S 

1  c 

c 
c 

LIMITATIONS  -  NOT  MORE  THAN  TWO  SWITCHINGS  AT  TIMES  TCI 

AND  TC2 

c 

r 

u 

INPUT/OUTPUT  LIST 

C. 

c 

VC1 

OUTPUT  VARIABLE  NO  1 

ANY 

OUTPUT  VAR 

c 

V02 

OUTPUT  VARIABLE  NO  2 

ANY 

OUTPUT 

VAR 

c 

MP\ 

INPUT  VARIABLE  NO  A1 

ANY 

INPUT 

VAR 

l  c 

VA2 

INPUT  VARIABLE  NO  A2 

ANY 

INPUT 

VAR 

\  c 

VB1 

INPUT  VARIABLE  NO  B1 

ANY 

INPUT 

VAR 

\  c 

VB2 

INPUT  VARIABLE  NO  B2 

ANY 

INPUT 

VAR 

\  c 

SMI 

SWITCN  CONTROL  INITIAL  VALUE 

— 

INPUT 

PARAM 

\  c 

*1.  VO-VB 

\  C 

»0.  VO-VA 

\  c 

tc: 

TIME  FOR  FIRST  SWITCH 

SECS 

INPUT 

PARAM 

\  c 

TC2 

TINE  FOR  SECOND  CWITCH 

SECS 

INPUT 

PARAM 

\J 

ITC2.GT.TC1) 

i 

• 

■n 

INE.GT.TCl.  AND.  TIME. LT.TC2) 


SW-ABS  ISHl-l.) 


vowai 

VC2=VA2 


Sw.GT.0.5) 


VOl-VBl 


(SW.GT.O.S1 


V02=VB2 


Table  96:  FLOWCHART  FOR  SUBROUTINE  SY 


\  SUWOUTIKE  STIV8t.V82.V93.V«l.Vfl8.Vft3.VJl.VK.V»3,Sm.TCl.Tca 

- — — , - - - 

_ I _ 

PURPOSE  -  T9  PROVIDE  a  SWITCH  CflWflNEKT  FAN  DftfE  VARIABLES 


NCTHflU  -  SEE  COSING 


WITTEN  IT  -  afiflK  LUSTS 


LATEST  REVISION  NOV  75 


lihitations  -  not  hk  than  two  switchings  aT  tikes  tci  aw  Tea 


input/output  list 


V81 

OUTPUT  VARIAILE  N8  1 

ANT 

OUTPUT  VAR 

V02 

OUTPUT  VARIABLE  WJ  2 

ANT 

OUTPUT  VAR 

VOS 

OUTPUT  VARIABLE  NO  3 

ANT 

OUTPUT  VAR 

vai 

INPUT  VARIAM  NO  Al 

ANT 

ItTUT  VAR 

vaz 

INPUT  VARIABLE  NO  A2 

ANT 

INPUT  VAR 

V«3 

INPUT  VARIABLE  NO  AS 

ANT 

INPUT  VAR 

VBl 

INPUT  VARIABLE  NO  B1 

ANT 

INPUT  VAR 

VB2 

INPUT  VARIAM  NO  B2 

ANT 

INPUT  VAR 

vbs 

tirur  variaSe  no  bs 

ANT 

INPUT  VAR 

SHI 

SWITCH  CONTROL  INITIAL  VAUUE 

— 

INPUT  PflRAN 

TINE  FOB  FIRST  SWITCH 
TINE  FAR  5ECON#  SWITCH 
tTC2.GT.TCl) 


input  panm 

IWRJT  P(W»t 


SU»5W1 

vn«vai 

vsa>va2 

vB3»v«3 


TWE-DT.TCl.flM.TlHE.LT.TC 


S»*AB3  tSWl-1.) 


IF 

N.GT.O.SL, 


Table  97:  FLOWCHART  FOR  SUBROUTINE  SZ 


\  SUIROUTINE  SI  (VS1 , V02.  V63.  YIN. Vfit .  YS2. YR3.  vm.  VS1 .  YC. VS3. YM.  / 
\i  3wt.TCt.Tca _ _ _ / 


w 

<- 

purpose  - 

to  mevicE  a  switch  cs^oneht  fsn  four  varmules 

c 

c 

c 

fd>3J  - 

see  coni ns 

c 

c 

t 

wimt*  pr  -  am*  turra 

LATEST  REVISION 

NOV  73 

t 

c 

r 

UMTHTIONS  -  NOT  HOSE  THAN  THO  SHITCHINC3  ST  TINES  TCl 

AN3  TC2 

c 

input/output  ust 

c 

vm 

output  v«ifl»f  no  i 
output  variaO  no  a 

AMY 

OUTPUT 

VAR 

c 

V32 

AMT 

OUTPUT 

VAR 

c 

VOS 

output  varia^x  no  3 

AMT 

OUTPUT 

VAR 

c 

von 

OUTPUT  VARIP1LE  NO  n 

ANT 

OUTPUT 

VAR 

r 

vm 

INPUT  VflRIAN  «  NO  A1 

mvr  vflRtfiix  no  re 

I*VT  YflNi.UU  NO  S3 

PKX 

IfTUT 

YflR 

VR2 

AMT 

iwnn 

VAT 

c 

vra 

pxr 

IHPUT 

VAR 

r 

vm 

IWUT  VARIANT  NO  m 

AMT 

ItPUT 

VAR 

c 

VII 

INVT  YARIA**  NO  II 

AMT 

INPUT 

VAR 

c 

VC 

INPUT  VflWfiBi  NO  C 

input  vmiftal  no  s? 

ANT 

iWUT 

VAR 

c 

vn 

ANT 

iirur 

VAR 

„  c 

VM 

DPUT  VflWfili  NO  M 

ANT 

INPUT 

VAR 

t 

/ 


f 

y 


'  C 

c 

C 

SMI 

sunch  CONTROL  INITIft.  Wilt 
-t.  V8-V* 

-0.  V5-VA 

_ 

INPUT 

FflRW  ^ 

v- 

TCI 

TI*  for  first  switch 

SECS 

INRJT 

PAR* 

V 

.  c 

TC2 

TINE  rm  SOCDHB  SNITCH 
ITC2.GT.TCH 

SECS 

INPUT 

PAR* 

SMHll 


voi-vai 

voa-vsj 

V83-V33 

vw»vm 


n 
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Table  V:  FLOWCHART  FOR  SUBROUTINE  SZ  (CONCLUDED) 


I  WO-VB 


Table  93:  FLOWCHART  FOR  SUBROUTINE  S2 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


I 


OF  3 

OUTPUTS  « 


H^'Sn'oS  ~  SUh  OF  For CES 

S££g£g-r  iSftf  ^^i5- 

L«*D1.F11  -  FORCES  I(*>VfT  port  , 

vpgft®1  ^ 2 


■x  ,  »  1  ‘  *  r  *  1 

flu.flm.ANi 

FJC2.FT2.FI2 

HL2, 

ST: 


SUM  FORCES 


HRY  1S78 


Table  100: 


FLOWCHART  FOR  SUBROUTINE  TA 


308 


nnoonnnnnnriri  onoooo 


Table  101:  FLOWCHART  FOR  SUBROUTINE  TB 


\  SUBROUTINE  TB«2T.B2T.fi2,Ba  / 


PU*,83€  "  I?  ?r£IDE  CfFflBILITY  TO  INPUT  VARIABLES  AS  FUNCTIONS 
Or  Tint 

METHOD  -  TABLE  LOOK-UP  USES  i  DEGREE  INTERPOLATION 
AND  WES  NOT  OSHO  JXTFW^ieK 


WRITTEN  BT  -  ADAH  LLOTD 


LATEST  REVISION  KSV  75 


IJFUT/OUTPUT  LIST 


aar 

B2T 


PS 

B2 


TABULAR  INPUT  OF  VARIABLE  PS  AKT 

JggjLflR  IfWT  OF  VAHIAOE  B2  MT 

ABOVE  TABLES  ARE  ALL  CHE  DIKENSIOWL 
WITH  TIME  AS  INDEPENDENT  VARIAME 
VARIABLE  A2  ANT 

VARIABLE  .32  ANT 


INPUT  TABLE 
IfttJT  TA®_£ 


OUTPUT  VAR 
OUTPUT  VAR 


T 


NA-A2TE2J  ~| 

NB-B2TC3 

^■TJUIl  ITIME.A2T  PC  .A2T 0»*C  ,  l.-NA) 


K-TSLUi  (TIME.B2TC10 
RETURN 


.32TWIMC  ;  I.  -NJD 


3C9 


Table  102:  FLOWCHART  FOR  SUBROUTINE  TBLl 
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Table  104:  FLOWCHART  FOR  SUBROUTINE  TBLU1 


c 
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\naCT19H  TH.U1  Qq.X.rv.HW.tft) 


atrcsE 


tMLffgfWM.Taiu  semen  ws  msnunoifw  ■’uuNflMa 
«?  uMhm  mm  m  \  rKotKNPtM 

s£T 


Vl. 

USS6£ 

BlNOfSIS*  XUTO.F1  WO 

5  l&Wi  i&T 
-g :  SBW&SF. 

HX  ** 


•  rayjMP&r rim  bwab  wr  ■ 


«umrr  rnwef 


■  wau  ar  rsur  stma*  mb  tMTxnMi<uiiM 

“S8S  ? :  IfflMW  '  ' 


WCNWWT  DISH 


(SBWCNOtj.x.ii&.ta,  n  .pj.ra 


Tab-e  104:  FLOWCHART  FOR  SUBROUTINE  T8LU1  (CONCLUDED) 


vL—  ~jWDWiArr> 

jgLM  »  TB»uxi,x.n.wnT7r 


am 


TabU»  105:  FLOWCHART  FOR  SUBROUTINE  TBLU2  (CONTINUED) 


Table  105:  FLOWCHART  FOR  SUBROUTINE  TBLU2  (CONCLUDED) 


Table  107:  FLOWCHART  FOR  SUBROUTINE  TD 


sunwutine  T3tv.v«.ry,p.pi,D*.n,M,r5,(W.l!Wj,iW.,  ~7 

\i  T<w.TflM».rfflH.r».pi8r.Mffrvm.ra.Ti.xxl.zii.xn.^rn/ 


E 


VERSION  l. 

“  * 


_  APRIL  77.  1877 

sr  *reej»  nicii  sarr  eouations  or  notion  i/fr 


mmmgufnflS 

linear  velocities  —  Msr  axes 

P.PE.IP  -  X  flilS  ANGULAR  vgjJCrTT. 

R.RI  "  “ 

■ft! 


ACCa.iKT  CONTROL  FT/SEC 


,INT  CONTROL. 
I  NT  CONTROL. 


R.R1.IR  -  2  fl(IS  ANGULAR 
■  ANGLES  —  EARTH  TO  tOt 

.R0LS.IR0L  -  RT ‘  - 

,rmjfiu,rrm  -  ru . 
mulm  —  earth  axes 

TS  -  T  AXIS  LINEAR  VELOCITT.  FT/SEC 

ANGULAR  ACCELERATION  -  WIT  AXES 

PHOT  -  X  AXIS  ANGULAR  ACSJJMTION.  DES/SEC2 

ROOT  -  Z  AXIS  ANGULAR  fESBWTION.  5eS/5€3 

mmmm  1(#VTJ  mm. 

UNEAR  ACCELERATION  —  MIT  AXES 

VBBT  -  T  AXIS  UNEAR  ACCELERATION.  FT/SEC2 


■flft 


:i&#or 


SUJG-fTS 


SLUG-FI2 


a 

V  c 


HRITTEN  IT 


—  EARTH  TO  WIT  -  PI7D1 
-  Pnot  ANGLE.  AEG 
H.K.  MAU 


HANOI  1877 


CPaC8SP!T<«PB 

5P-SIA  tpn«RPB 

OW35IWL^PB 

SRtSMMUtiPB 
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Table  107:  FLOWCHART  FOR  SUBROUTINE  TO  (CONTINUED) 


j 


M-rn/nn«m 


LiOla  108:  FLOWCHART  FOR  SUBROUTINE  TERRA 


f 


VWiM 

mnts 


i 


DISCRETE 


r 


FUNCTION  TERRA  (X,AM0DE,ANR,D,M,Z 


Hi  •  AMR 
N2  •  ANR*3 


-»(  RETURN 


IMS  *  XE/D 
1  -  IM1+4 


TERRA  -  Z{4)*N 


TERRA  *  Z(M2)«M 


XS  •  XHI*0 
BE  •  XE-XZ 

TERRA  •  (UU*1W(I))WD*I(I))*H 


J. 


-|  TERRAKI 


I  TERRA  »  .S*N*(l.-COS(6,2832*XE/D)) 


nnnnnnnnn'  irioooooooooocmoon 


Table  109:  FLOWCHART  FOR  SUBROUTINE  TF 


\  subroutine  tf txi.'xi3gT.ixi.F5.fsagr.if3.r:s.zo.-i.°o.?ii 

_ i _ 

PURPOSE  -  T0  SIMULATE  P  SECOND  5TDER  TRANSFER  FUNCTION  *ITH 
FIRST  ORDER  KjMERRTOR 


FO 

FIN 


METHOD  -  SELF  EXPLANATCRT 


£i«S  -  £0 


S  -  ?1-S»P0 


LIMITATIONS  -  NONE 

WITTEN  3T  AGAN  LL?v-3  LATEST  REVISION  «JV  7S 


INPUT/ OUTPUT  LIST 


XI 

INTERMEDIATE  STATE  VARIABLE 

ANT 

X'TPUT  STAT 

X1B07 

STATE  VfRIABL £  DERIVATIVE 

ANT 

OUTPUT  STAT 

1X1 

INTEGRATOR  CONTROL 

— 

PROGRAM  VAR 

FO 

TRANSFER  FUNCTION  STRUT 

ANT 

OUTPUT  STAT' 

FODOT 

IFO 

TRANSFER  FUNCTION  OUTPUT  DERIV. 
integrator  CONTROL 

ANT 

OUTPUT  STAT 
PROGRAM  'iAR 

FIN  TRANSFER  FUNCTION  INPUT  ANT  INPUT  VAR 

~  *  .  -i  -  - 


i _ 

r 

10 

NUMERATOR  COEFF  ICIEST 

ANT 

IIPUT 

VAR 

\ 

r 

21 

numerator  coefficient 

ANT 

IWJT 

VAR 

•v 

PO 

DENOMINATOR  COEFFICIENT 

1/SEC2 

IWJT 

VAR 

; 

r 

PI 

DENOMINATOR  COEFFICIENT 

l/S££ 

INPUT 

VAR 

UJU.)  UJUJ 


nnnnnnnnnnnn 


Table  110:  FLOWCHART  FOR  SUBROUTINE  TG 


SUBROUTINE  TG  lFX.FT.Ft.TX.TT.TZ.TH.  GflMX.GfiMT.GfiMI.X0.T0.20) 


VERSION  1.  REVISED.  SEPT  20  1S766 

PURPOSE:  TRANSFORM  EHGIKE  7WTUST  INTO  3CBT  AXIS  FORCES  AND 

TOROUES 

CALL  StglEHCE:  . 

FX.FY.Fz'  '  -xXz  AXIS  FORCES  ABOUT  C.G..  LBS 
TX.TY.TZ  -X.T.Z  AXIS  TORSUES.  FT-LBS 

TK  -ENGINE  TWJST  IN  LBS 

Gfftt.GAHr.GflMZ  -X.T.Z  AXIS  DIRECTION  COSINES 
XQ.TO.ZO  -  TFflUST  LOCATION  OPPONENTS  .  FT 
DESIGNED  BT  A.M.  HRRREN  SEPT.  1376 


FX»  TH-GTX 
FT^  TH«&*-.r 
FT*  TA«G  *^ 

TX»  YO-F!  -  ZOFY 
TY»  ZU  d  X  -  XG-FZ 
tz=>  xo-fr  -  ro-fx 
return 


Table  111:  FLOWCHART  FOR  SUBROUTINE  TK  (CONTINUED) 


rf*  8a  3  3  3 


Table  112:  FLOWCHART  FOR  SUBROUTINE  TL  (CONCLUDED) 


CT<3I(YrtW*» 

sr-siNnrawvB 
x>cr<r><^  or«Sit*a»aPOT  «*< 


Table  113:  FLOWCHART  FOR  SUBROUTINE  TR 


\  SUBROUTINE  TRPXE.PTE.PZE.PXB.PTB.PZB.ROL.PIT.YAW  / 
VERSION  l  JO-E3  1377 

TRANSFORM  VECTOR  QUANTITIES  FROM  BOUT  /WES  TO  EARTH  AXES 
MATRIX  MULTIPLICATION 


C — PURPOSE 
C 

C— METHOD 
C 

C— CALL  SEQUENCE 

C  . in  OUTPUTS 

C  PXE.PTE.PZE 

C  INPUTS 

C  PXB.PYB.PZB 

C  ROL.PIT.YAM 

C - WRITTEN  BT  J.J.HCAVOT 

C  CONVERT  FROM  DEGREES  TO  RADIANS 


VECTOR  QUANTITIES  ALONG  EARTH  X.Y  AND  Z  AXES 


VECTOR  COMPONENS  OF  BODY'  COORDINATE  SYSTEM 
BODY  ROLL. PITCH. AND  YAH  ANGLES 


R3L-R0L/57.3 

PIT-PIT/S7.3 

YAH-YAN/57.3 

_L 


PXE-PXB-CCS  PIT)  -COS  CYAN)  ♦fYB-  I-COS  tROU  -SIN  CYAN)  *SIN  (ROU  - 

2  SIN  (PIT)  -COS  OfftO )  *PZB"  (SIN  (YAW  -SIN  (ROU  -C&3  (YAH1  -COS  (ROD  « 

3  SIN  (PIT) ) 


(C— EARTH  LATERAL  COMPONENT) 
- 1 - 


FYE*-PXB»*COS  PIT)  -SIN  (TFN)  -PYB-  (COS  (TAW  «COS  WSJ  ♦SIN  (TAW  * 

2  SIN  (ROU  -SIN  PIT)  3  -PZ8«  (-COS  (TAW  -SIN  (ROU  *SIN  (YAW  -COS  (ROU  < 

3  SIN  PIT)) 


CE 


I 


-EARTH  VERTICAL  COMPONENT") 


PZ^PXB-SiNPIT)  -PYB-C9S  PIT)  -SIN  (ROU  -PZB-COS  P.OU  -COS  PIT)  1 


B30 


PAGE  i 


Table  114:  FLOWCHART  FOR  SUBROUTINE  TRUNK  (CONCLUDED) 


Table  115:  FLOWCHART  FOR  SUBROUTINE  TS 


_ , 

n  Win  H*U  l  UJ78 

wmn  ir  i.i.xlor 

hi w*  -  tb  Knur  m  amtc  ms  kls  dim  i 
c  «m»  -  as  vaui«  i 


COOTWTS  9mm  _ 

fiwo.3i  -a»  ctnsTwr  rm  am.  rr-u/u/wo* 

swaws 


wwHtt8»tam«'»i.:»sna»3.  1 


u 


— M  ao  5  M.m  I 

r 


s  •resn«TtsTt*»iH(n 


nt»<T>i(a«(yma*t.)*sTw  a  *s. 


5SP 


PlW  t 
TS 
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Table  11*5;  FLOWChART  FOR  SUBROUTINE  TS  (CONTINUED) 


Table  115:  FLOWCHART  FOR  SUBROUTINE  TS  (CONTINUED) 


Table  115:  FLOWCHART  FOR  SUBROUTINE  TS  (CONTINUED) 


IE  TRN<  ELOWT  IDOBIfl*  AM  fi W9  fW  F\I  Sfff 


MUi  i  t*l 


Table  115:  FLOWCHART  FOR  SUBROUTINE  TS  (CONCLUDED) 


uuououuuuuuuuuoouuuuuuuuuu 


Table  117:  FLOWCHART  FOR  SUBROUTINE  TZ 


\  SUBROUTINE  TZ  (Xl.XlDOT,  IXI.X2.X2D0T.  IX2,F0.FIN.Z0,Zl.Z2.P0,Pl)  7 

.  ~  I  :.::- 

PURPOSE  -  T5  SIMULATE  SECOND  ORDER  TRANSFER  FUNCTION  WITH  SECOND  ORD 
NUMERATOR  2 

FO  Z2«S  *  Zl*S  *  ZO 


FIN 


S  *  Pt«S  +  PO 


METHOD  -  SELF  EXPLANATORY 


LIMITATIONS  -  NONE 


WRITTEN  3Y  -  ADAH  LLOYD 


LATEST  REVISION  NOV  75 


INPUT/OUTPUT  LIST 

XI  FIRST  STATE  VARIABLE 

XI  DOT  FIRST  STATE  VARIABLE  DERIVATIVE 

1X1  INTEGRATOR  CONTROL 

X2  SECOND  STATE  VARIABLE 

X2D0T  SECOND  STATE  VARIABLE  DERIVATIVE 

1X2  INTEGRATOR  CONTROL 


ANY 

OUTPUT  STATE 

ANY 

OUTPUT  STATE 

PROGRAM  VAR 

ANY 

OUTPUT  STATE 

ANY 

OUTPUT  STATE 

PROGRAM  VAR 

ANY 

OUTPUT  VAR 

_ i _ 

C 

FIN 

TRANSFER  FUNCTION  INPUT 

INPUT 

VAR 

c 

ZO 

NUMERATOR  COEFFICIENT 

INPUT 

PARfiM 

c 

Z1 

NUMERATOR  COL"  ICS  ENT 

INPUT 

PARAM 

c 

Z2 

NUMERATOR  COEFFICIENT 

ANY 

INPUT 

PARAM 

c 

PO 

DENOMINATOR  COEFFICIENT 

1/SEC2 

INPUT 

PARAM 

C 

PI 

DENOMINATOR 

1/SEC 

INPUT 

PARAM  „ 

FO»X2*f IN«Z2 


Table  118:  FLOWCHART  FOR  SUBROUTINE  VA 


VERSUS*  l. 


umm*  VAU9. V9.HS.P9. 96.A0, 3iN9,M}.f949,CRS,  JAS.*..ALP. 
VMS,  JETA,  HP, IF,  QJ.CV.Q4.SiSZ.  SCSN.  AJAR.  »ACH.  U.  V,t4. 

p.  a.Kz.Pi-.  *L  bin.  v&.  als.  s.  vu.  vm.  m.  pm.  cm.  ahTbio 

rz 


amrTE  «Re  varuules  ns*  states 


HAT  It  |  87# 


airmrr  variaxxj.  iiryTs  -  - 

ue.vfl.we  -x.t.I  *esT  axis  linear  VELgcrrr  with  mixs-ft/sfc 

po.se.fse  -x.T.i  bbit  axis  angular  Sites  urm  wins.  jes/sec 

»nse  -iwicm»  f*  aejSISsces  m  n*@(ts  -  sir 

SNA.  R4fl  -9  AKI  IS  ANGULAR  RATE  OUSTS.  3ES/S E 

KK  vm'MLtt  -  -  _ _ 

-  - Ti«  MSIlCS  FOR  STMJUTT  «J3TRRNSFflm 

of  attaoHa  i esr  am  stawutt  axes,  ses 

'  -TT 


CAS, SAS 

al.alp 

VMS 


5!^ 

EU.EV.EW 


I1**  X  STAliLl'T  AXIS  velocities, ft/sec  UINQOIflN 
-!«lt  PQmiUATION  VELOCIT I OOMiHCWSi ONAU 
-X.T.T  40BT  AXIS  ACCELERATION  TOW  FAR  USen.VBOT. 

wiot  solutions.  ft/sec-s 

STANiWB  ATNOSPHERE  VARIABLES  -  - 

siez  -aLs  awsiTTfsAT 

seas 

SMS  _  _ 

XHACH  -4SRQ4  liBUCR 

STATE  VARUULES  ' 

u.v.w  -x.t.i  sen  axis  linear  velocities.  ft/see 


-owressixe  urwesic  pRCssme,  uvrr«»»2 

-8TMNIC  °RESSUA£  TIRES  AEFEROifcE  AREA.  JS 


P.S.R 

z  _  . 

PIT  RBL  -PITCH  AM  ABU 

IN1ICAT0R  FUNCTION  FflR  AERO  Ff 


I 


ai  atm 


1IR 


CMSTMT3  - 
VS 

f5 

WIKI  STATES 
UH.VM.H4 


-KET  AXIS.BINEKSION*. 

-BOUT  AXIS.  WNfllNEKSION*. 
-STAlrijTt  AXIS.  BlNEKSIOM. 
-STMIUTT  AXIS,  MMIREHSii 


KS/SEC 

....  I  TO'IMT  AXIS  EULER  ANGLES.  3C 
AM  AflHWTS 


AIMKES,  FT/XC 
—  -  OF  ATTROC  BEE 

—a 


X.T.J  SOfT  «I5  HIM  VELOCniES^  JT/gC 


BflBT  AXIS  WIM  ANGULAR 


PM  CM  AN  -X  X  t 

UUltRT*  FUNCTION  F«  8EGREE3  OF  FREES*  IEOH 
816  -  a  -TWO  JOF  LflNfilTVBIN*.  15,0 

3  -THREE  DDF  L9NEITU91NRL  U.U.GS 

4  -FOUR  SflF  LATERAL  tv. P,R>  *  LflNSJ.  8B 

5  -•MiejorLATSWL  (V.P.W 

_ S  -«IX  30F  FUU.  NOBEL  W.Y.W.P.8.FD 

MITTEN  ST  WWIH8BS  WftU  NAT  1877 


6/SEE 


lHmnjmim 


PAGE  1 
YR 
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Table  118:  FLOWCHART  FOR  SUBROUTINE  VA  (CONTINUED} 


Table  178:  FLOWCHART  FOR  SUBROUTINE  VA  (CONCLUDED) 


IF  ^ 
UN.C0.1.L. 


.  JF 

3IH.E8.3.) 


if*  ma  -  un /vs  i 


tp»  ay  -  vsi/vsl 


ISVSVS^ 
JfcffliSflR  aagRffll 


?  cj  1  ■  “  —  ■« 

'm  of  Ain  Bw»m  iktis  am  saws  velocity 


«-  2/3000.  *t.  - -  - 

if *_HIH0  ««0  (t .  IZ) .  3U> 

Qtm-  OPSi  s 
M«CH«  YWVflJ 
X*f2-  /JVSX—S 

_^9>  DPS*!!.*  XXg».a«ll.»  XH?«.C233)) 


‘'CfxmJr. 


OCffy  0PS«U.93a  -.772/»«a— 3  *  .flisS/XWfrj^T 


t  4 1 


SC5E3I 


Table  120:  FLOWCHART  FOR  SUBROUTINE  VLX 


355 


356 


Table  121:  FLOWCHART  FOR  SUBROUTINE  VPRINB 


357 


flowchart  for  subroutine 


V?RINB  (CONCLUDED) 


j*ms(IXtnr.aoT 


Table  122:  FLOWCHART  FOR  SUBROUTINE  VPRINT  (CONCLUDED) 


Table  123:  FLOWCHART  FOR  SUBROUTINE  WS 


\  sumrTiNc  to rog.uHs.vws.we. tx.HRH.fltT.pn.ain 


X 


vgai*  2. 

PURPOSE  SII 
METHOD  -Mil 
NODI 


MARCH  31  1378 
irr-yma  rarmntn 


^ \ 


UDE-OF  SQ-EEET  IS 


_  _ L  _,u. _ TOR  TO  REFLECT  THE 

CHANCE  IN  HIND  HIIK.CTlTC5f.NWa  YCCTW  WSSJMEB 
PWm  TO  THE  I3W.1W  PUNg.THE  WND  MAONITeOE  MODIFIED 

IT  SHEW  IS  RESOLVED  HONG  THE  BIDKHT  COORDIWTTES  WORTH 

m»  east  lfwj  “side  and  nwHsfimj  into  bsdt  axes. 


SLPflg  — 

no  -ONE  HiNENStflNALTAJLE  OF  MINS  S«AR  FACTOR  AS  A 
-njHCTION  OF  AIRPLANE  CO  fl-TITUDE 

LINEAR  VoSciTIEs"— '«irr  SCES 
UH3  -X  AXIS  STEADT  OR  SMEAR  HIND  VELOCITT.  FT/SEC 

VMS  -T  AXIS  STEADT  OR  SHEAR  MIND  VELOCITY.  FT/SEC 

MtS  -I  AXIS  STOUT  OR  SHEAR  HIND  VELOCITT.  FT/SEC 

INPUTS  n— 

-HIND  MAGNITUDE  AT  SO  FEET  CT8MERJ .  FT/SEC 
-ANGLE  SETHEER  THE  HIND  VECTOR  AND  RUMAT  CENTERLINE,  DEB 
-AIRPLANE  CO  ALTITUDE.  FT 

-PITCH  ANCLE  EARTH  TO  MT.  DEB 
INDICATOR  FOR  STEADT  OR  SHEAR  HIND 

MU*  O  SHEAR  Hi NO, TAILE  LOOKUP  FACT*  y 

»1  STEADY  HIND.  FWCTOR»l. _ / 


IK 


H.T 

PIT 


r 


WRITTEN  IT  WINDER  UAHI 
SET DEFAULTS 


MAT  1377 


0 


Table  124:  FLOW  FOR  SUBROUTINE  XP 


SUBROUTINE  XPrT,P2,02.R2,Pl,Qi.Rl) 


VERSION:  1.  REVISED:  JUNE  10  1976 

PURPOSE:  PERFORM  STATIC  TRANSFORMATION  ON  TWEE  VECTOR 

CALL  SEQUENCE-. 

T  -  TRANSFORATION  MATRIX 

P2.Q2.R2  -  OUTPUT  VECTOR  COMPONENTS 

Pl.Qi.m  -  INPUT  VECTOR  COMPONENTS 

DESIGNED  BY;  J.D.  BURROUGHS  JUNE  1376 


P2-=pi»*t  am  *oi«t  a  n  »ri«t  u2) 
Q2=Pl«T  (13)  -HU-T  (HI)  +RI-T  US) 
R2-=Pl«T  (16)  *01«T  (17)  *Rl«T  (18) 
RETURN 


PAGE  t 
XP 


MW*****' 


nnnnnnn 
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FLOWCHART  FOR  SUBROUTINE  YC  (CONCLUDED) 


r 


i 

I 

369  pqr,r  ?  j 

I 
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